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Abstract. The question of self-adjoint realizations of sign-indefinite second 
order differential operators is discussed in terms of a model problem. Operators 
of the type — sgn(x) are generalized to finite, not necessarily compact, 
metric graphs. All self-adjoint realizations are parametrized using methods 
from extension theory. The spectral and scattering theory of the self-adjoint 
realizations are studied in detail. 



1. Introduction 
Differential expressions of the form 

(1) - divA(-) grad with A(x) = < 



-1, xeSl + , 
-1, xeSl_, 



for domains = £1 + U £!_ C H d appear in different contexts. The Poisson problem 
related to this operator is of physical relevance. It appears in the mathematical 
description of light propagation through regions with thresholds of materials with 
negative and positive refraction indices, see [5] and the references therein. The 
quasi-static limit of the Maxwell equation yields the mentioned Poisson problem. 
Another field of application is solid-state physics. In the effective mass approxi- 
mation the effective mass tensor can be negative as well, see for example [T] and 
the references therein. This yields indefinite differential operators, which involve 
expressions of the form ([1} and in this context the Schrodinger equation is consid- 
ered. 

An approach that has turned out to be very fruitful for the study of differential 
equations which involve the expression — div A(-) grad with strictly positive coeffi- 
cient matrices A(-), is to consider self-adjoint operators associated with these dif- 
ferential expression and their spectral resolution. The question is, if it is possible to 
find self-adjoint realizations also for sign-indefinite coefficients A(-). Unbounded 
operators with sign-changing coefficients in the highest order terms are not well 
studied. To make a starting point a new model problem is introduced: Laplace 
operators on finite metric graphs with changing sign. The aim of this work is to 
present generalizations of the operator 

(2) ^-sgn(x)-^- with sgn(a;) = 

ax ax 

to finite metric graphs. This initial problem can be tackled using extension the- 
ory. This approach allows straight forward generalizations from intervals to finite 




2010 Mathematics Subject Classification. Primary 34B45, Secondary 47B25, 34L05, 35P20, 
35P25, 81U15. 

Key words and phrases. Indefinite differential operators on metric graphs, self-adjointness, 
spectral theory, scattering theory. 

1 



2 



AMRU HUSSEIN 



metric graphs. It yields a formal operator which admits the block operator matrix 
representation 

_ [-A " 
T ~ [ +Aj ' 

where —A denotes the one dimensional Laplace operator, and all self-adjoint real- 
isations are completely characterized. On finite metric graphs this gives a model, 
which on the one hand is still explicitly solvable, but on the other hand allows to 
describe more complicated geometries than the one dimensional case does, which 
deals only with intervals. A similar procedure has been performed already very 
successfully for the sign-definite one dimensional Laplacian. This has become fa- 
mous under the name "quantum graphs". For further references on this topic see for 
example [3J and the references given therein. It turns out that there is a one-to-one 
correspondence between the self-adjoint realizations of the Laplacian on graphs and 
the self-adjoint realizations for the model-operator discussed here. 

The aim is to construct self-adjoint realizations of expressions of the form (J]} 
and ^ in terms of extension theory and to develop a better understanding of the 
spectral and scattering properties of sign-indefinite differential operators. Some 
aspects of the problem are included in the general setting of Sturm-Liouville theory, 
but they have been considered to have no practical relevance, as remarked in |29L 
Section 17.E]. 

The motivation to consider expressions of the form (JTJ in the present paper 
is caused mainly by the fast development in the field of so called metamateri- 
als - materials with negative optical refraction index. This has been initiated by 
V. G. Veselago, see [25]) who was the first to consider - hypothetically - the effect 
of a negative refraction index. He had predicted that this feature would yield new 
effects including for example backward waves, an inverse law of refraction and an 
inverse Doppler effect. Today metamaterials can be constructed as artificial ma- 
terials, whose optical properties are made useful for applications. An overview of 
recent developments is given in |27| . recall also the article [5] and the references 
given therein. 

Note that the model problem ((2]) is part of a range of problems related to the 
ordinary differential expressions 

d , . d 

and — sgn(x)— sgn(x) — . 

dx dx 

Each is exemplary for certain difficulties and methods used. Usually the operators 
written in the first line are considered in Hilbert spaces, whereas the operators in 
the second line are considered in Krein spaces, see for example [6, 32J. The operator 
— sgn(x)-^ has applications in the effective mass approximation too, compare [33]. 

This note is organized as follows: first the basic concepts are introduced and the 
notation is fixed. The subsequent section is devoted to the characterization of all 
self-adjoint boundary conditions. In short the corresponding question for Krein- 
space-self-adjointness is discussed as well. Section [5] puts the question of self- 
adjoint extensions into the context of indefinite quadratic forms, and in Section [BJ 
the operators discussed are put into the general context of extension theory. Section 
Ogives explicit formulae for eigenvalues, resonances and resolvents. In Section[5]the 
scattering properties of the system are discussed. The wave operators are computed 
as well as the scattering matrix. The scattering matrix can be computed in certain 
cases in terms of a generalized star product. 



dx 2 ' 
d 2 
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The content of this note has developed in parallel to the work on the article 
[TB] which is in preparation, and it is part of the author's PhD thesis, see [T51 
Chapter 4]. Indefinite operators on graphs provide an illustration and can serve 
as a source of examples for indefinite operators of the form ([T]) on manifolds and 
intervals involving more general coefficients A(-). 

Acknowledgements. I would like to thank David Krejcink for communicating the 
problem of sign-indefinite differential operators, its physical background, helpful 
discussions as well as remarks and especially for kind hospitality during my visit to 
the Doppler Institute in Prague in February 2012. I am grateful as well to Vadim 
Kostrykin and to Stephan Schmitz from whom I benefited in many ways during the 
still ongoing work on the article [J2| . 

2. Basic concepts 

Instead of considering the operator given in @ in the Hilbert space L 2 (R, dx), 
one considers in the equivalent space 

L 2 ([0,oo),d Xl ) L 2 ([0,oo),dx 2 ) (= L 2 (R,dx)) 

the formal differential operator 

TU = (+<,-<)• 

This allows to define a closed symmetric operator with equal deficiency indices 
(2,2). Hence there exist self-adjoint extensions. The operator r is generalized 
from intervals to finite metric graphs which delivers indefinite symmetric second 
order differential operators with finite deficiency indices. In the following the ba- 
sic concepts are presented. The notation is mainly borrowed from the works of 
V. Kostrykin and R. Schrader on Laplacians on metric graphs, see for example |19| . 

2.1. Finite metric graphs. Here and in the following a graph is a 4-tuple 

Q = (V,l,£,d), 

where V denotes the set of vertices, I the set of internal edges and £ the set 
of external edges, where the set £ U T is summed up in the notion edges. The 
boundary map d assigns to each internal edge i G I an ordered pair of vertices 
d(i) = (wi,«2) G V x V, where v\ is called its initial vertex and V2 its terminal 
vertex. Each external edge e G £ is mapped by d onto a single, its initial, vertex. 
A graph is called finite if |V| + |I| + \£ | < oo and a finite graph is called compact 
if £ = holds. 

The graph Q is endowed with the following metric structure. Each internal 
edge i e I is associated with an interval [0, a^] with > 0, such that its initial 
vertex corresponds to and its terminal vertex to at . Each external edge e G £ is 
associated to the half line [0,oo), such that 9(e) corresponds to 0. The numbers 
cti are called lengths of the internal edges i €l and they are summed up into the 

\T.\ 

vector a = {a^gx G R + . The 2-tuple (Q, a) consisting of a finite graph endowed 
with a metric structure defines a metric space called metric graph. The metric on 
(<7,a) is defined via minimal path lengths on it. 
One distinguishes two types of edges 

£ = £ + Q£- and 2 = 1 + Ul- , 

where £ + {£-) and T + are called the positive (negative) external edges and 

positive (negative) internal edges, respectively. The set £+UT + (£_ Ul_) denotes 
the positive ( negative ) edges. This partition defines two sub-graphs 

G-,G+CG, where G± = (V±,l±,£±,d \ {£±ux±) ) 
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with V± = d(£±Ul±). Let be a ± = {di} ie z ± , then (G±,&±) are metric graphs. 
For brevity set 

n = 2|I + | + and m = 2|I_| + \£-\. 

2.2. Function spaces. Given a finite metric graph (Q, a) one considers the Hilbert 
space 

H = H(£,l,a) = n e ®Hi, H £ = 0"H e , Hz = 0H l7 
where "Hj = L 2 (Ij) with 

[0,aj], ifjeX, 
JO, oo), iijeS. 
Again one has a decomposition into two subspaces 

H(£,I,&) = H(£+,l+,a i+ ) 8 %{£-,!-,&_). 
Any function (C/,a) — > C can be written as 

tp(xj) — i>j(x), where ipj : Ij — > C 



with 



One defines 



[0,aj], ifjeX, 
[0,oo), ifje£. 



p&i />00 

ip ■= y~] I ip(xi)dxi+^2 ip(x e )dx e , 
tex^ ees Jo 
where dxi and dx e refers to integration with respect to the Lebesgue measure on 
the intervals [0,aj] and [0, oo), respectively. 

By Wj with j e £ U I one denotes the set of all ipj € "Hj such that its derivative 
V>j is absolutely continuous and is square integrable. On the whole metric graph 
one considers the direct sum 

W = W(£,1)= W,. 

je£ ui 

By 2?j with j £ £UI denote the set of all ^ € "H j such that ^ and its derivative 
ip'j are absolutely continuous and tp" is square integrable. Let T>j denote the set of 
all elements ipj G T>j with 

V,(o) = o, V'(0) = 0, forje£, 

■0^(0) = 0, V'(0) = 0, ipj(a,j) = 0, ip'(a,j) = 0, for jel. 
In the space % one defines the subspaces 

V= V h and P°= Uj. 

jesux jesux 

Note that these spaces clearly decouple the edges of the graph. The coupling is 
going to be implemented in terms of boundary conditions. With the scalar products 
(•, -)w defined by 

((p,ip) w = (<p,i>)H + (<p',iP')h 

and (-, -) v given by 

M)<d = M)n + y,V>')« + (<fi",ip")n 
the spaces W and T>, respectively become themselves Hilbert spaces. 
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2.3. Differential operator, in this article the formal differential operator 
(3) {r^) j (x) 



d 2 



is considered. That is r acts as — A on (Q + ,a + ) and as +A on (C/_,a_), where 

d 2 



(4) 



(A^) j (a) 



j e £ ill, x e 



To emphasise the importance of the domains and to make the choice of the do- 
mains transparent, the notation is used to write a differential operator as 2-tuple, 
consisting of a differential expression and a domain. The operator 

T mi " = (T,2?o) 

is closed, densely defined and symmetric. Its adjoint in the Hilbert space % is the 
operator T max = (T min )*, 

T max = 

Since the operator T mm has equal deficiency indices 

, d-) = (d,d), where d = n + m, 
there are self-adjoint realizations T — T* of r with 

and all these realizations can be parametrized using boundary conditions. 

2.4. Space of boundary values. On the lines of [T^] one introduces the auxiliary 
Hilbert spaces 



JC± = 1C(£±,2±) = lCs i 



with JC £± = Cl £ ±l and JC^J T ±K Eventually one defines 

/C = /C+©/C_ and K 2 = K,@JC. 
For ip € T> one defines the vectors of boundary values 



{^(o)W + 

{^(0)}, e z + 
{^i(ai)}iex + . 

{^•(0)b eC+ 
{^•(0)} i6 x + 
{-^K)bex + , 



"{^•(o)W_' 

{^•(0)} ie z_ 
.{^j(oj)}jex_. 

{^•(0)be£_ 
{^(0)bex_ 



if 



The vector [tp] = tp © -0' is an element of K, = K © /C. 



3. Self-adjoint realizations 

To measure how far the maximal operator 7 imax is away from being self-adjoint 
one evaluates the form }(g +l g_) defined by 

kg + ,g-)U>,<P) ■= (T max ^> - faT™*?), for i>^eV, 

compare for example [141 Section 3]. Note that this form can be represented by a 
matrix in the finite dimensional Hilbert space K. , that is 

>(<?+,!?-) W></?) = ([#%+,S_)M>JC 2 , 
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where the matrix 



3(s+,a_) 






l n 





-1, 







is written with respect to the decomposition of K 2 = /C+ © JC- © /C+ © /C_; l n 
denotes the identity operator in the n-dimensional space /C + and l m the identity 
in the m-dimensional space /C_, respectively. By an abuse of notation one denotes 
by j(g,,g_) also the sesquilinear form defined by 3(g + .g_) on JC . Notice that 



-1a 



-•J(0+,s_) and ^(e+,e_) 

holds and therefore is a Hermitian symplectic matrix in /C 2 . 

of Z(g + ,g_) to the standard Hermitian symplectic matrix 



The relation 






In 




"In 







is provided by the similarity relation 



(5) 3g = Hn,mZ{g + ,g_) #n,m, where H n>r 





-1, 








1„ 





A subspace Ai C /C 2 is called maximal isotropic with respect to j/g +j g_) if 

i(e + ,e_)(6>6) = (£i J 3g + ,a_ 6) = for all ^^eM 

and if there are no proper subspaces M. C M. having this property. Any subspace 
M C K 2 with d\mM > d, can be parametrized as follows. Let A and B be linear 
maps in 1C. By (^4, B) one denotes the linear map from JC 2 = K © K to /C defined 
by (A, -B)(?7i © 772) = At]\ + where 771, 772 € K. One sets 

M = M(A,B), where M(A -B) := Ker(A, B). 

From ([5]) one deduces a one-to-one correspondence between subspaces in JC 2 which 
are maximal isotropic with respect to the standard Hermitian symplectic form, 
which is well studied, and subspaces in K 2 that are maximal isotropic with respect 
to the Hermitian symplectic form )(g +: g_) occurring here. 

Theorem 3.1. All self-adjoint extensions o/T mm are uniquely determined by sub- 
spaces Ai C fC 2 , which are maximal isotropic with respect to the Hermitian symplec- 
tic form )ig + ,g_) ■ All such maximal isotropic subspaces are given by M. = Ai(A, B), 
where A and B are linear maps in K,, which satisfy the two conditions 

(1) Rank(A, B) = n + m, that is the rank is maximal and 

(2) BJ n ^ m A* = AJ n , m B* , where J n ^ m : "" " 







-1, 



defines a map in 



/C = /C+©/C_. 

Any self-adjoint extension of T mln is given by T{A 1 B) = (r, Dom(T(A, B))) with 

Dom(T(A, B)) = {<i/> G V \ Af_ + B^_ = 0} , 
where A, B satisfy both conditions (1) and (2). 
Remark 3.2. 

(1) The condition Ai\) + Btp' — is equivalent to [tp] £ M(A, B). 
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(2) The operator —T(A, B) is the operator with positive and negative edges 
interchanged, but the coupling is implemented by the same boundary condi- 
tions at the vertices. 

(3) For — the operator T(A, B) is the self-adjoint operator plus or minus 
Laplace =fA(A, B) on Q±, where A(A,B) denotes the realization of A with 
domain Dom(A(A, B)) — {ip £ T> \ Atp + Bip' = 0}, compare for example 

(4) Observe that the second condition in Theorem can be re-formulated 
in terms of relations in finite dimensional Krein spaces. Consider the 
space K, equipped with the indefinite inner product [•, •] = (•, J n ,m')K- Then 



rt 4* 



BA^*>, where •<■*' denotes the Krein space adjoint. 



Proof of Theorem \3.1\ It is well known that all self-adjoint extensions of T mm , 
which is a closed symmetric operator with finite and equal deficency indices, is 
determined by subspaces A4 of /C 2 , which are maximal isotropic with respect to the 
Hermitian symplectic form see for example P3] and the references therein. 

To prove the parametrization let now M. be a maximal isotropic subspace of K 2 
with respect to the Hermitian symplectic form defined by Z(g + ,g_)- Then by §5§ the 
space H n m M. is a maximal isotropic subspace of K. 2 with respect to the standard 
Hermitian symplectic form defined by Zg- The space B. n ^ m M. can be parametrized 
by matrices Ao,Bq, see for example [19) . which satisfy the two conditions 

(1) Rank(Ao, -Bo) = n + m and 

(2) A B% = B A* . 

Hence the original space M. can be represented as 

"1„ 
-1, 



M = H n , m Ker (A , B ) = Ker ^A , B 
The matrices 

A = A and B = B a 
satisfy the two conditions formulated in Theorem 13. II 



In 

-1, 



□ 



Remark 3.3. The proof of Theorem VS. 1\ shows that if one has a self-adjoint Laplace 
operator — A(Aq,Bq) parametrized by matrices Aq and Bq, then T(A,B) is self- 
adjoint with A = Aq and B = B§J n ^ m . Ln particular there is a unique parametriza- 
tion in terms of unitary operators in K,. For Laplace operators it is known that for 
self-adjoint boundary conditions defined by operators Aq and Bq in K, there exists 
a unique unitary map U in K, such that equivalent boundary conditions are defined 
by 



1 



1 



A' = --(U-t) and B' Q = -(U + t).. 

see for example |14l Section 3]. Hence for T(A,B) self-adjoint one can give an 
equivalent parametrization in terms of the same unitary matrix U with 

A' = -l(U-t) and B' = ~(C/" + l) J n ,m- 

The parametrization by matrices A and B in Theorem 13.11 is not unique, since 
two operators T(A',B') and T(A,B) are equal if and only if the corresponding 
maximal isotropic subspaces Ai(A' , B') and A4(A, B) agree. Using the result from 
[231 Theorem 6], for self-adjoint Laplacians on finite metric graphs one obtains a 
unique parametrization. 
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Corollary 3.4. Let M.(A,B) C /C be maximal isotropic with respect to the Her- 
mitian symplectic form )(g + .g_). Then there exists a unique orthogonal projection 
P and a unique Hermitian operator L with P^LP^ — L, where P = 1 — P, such 
that with 

[i» 



A' = L + P, 



B' = P 







-1, 



one has M(A,B) = M(A' , B'), where P is the orthogonal projector onto 
Ker BJ n m C K,. 

Example 3.5. Let Q be the star graph consisting of two external edges, £+ = {1} 
and £ _ = {2}, glued together at one single vertex 9(1) = d(2), and consider 



A 



-1 





and B 





-1 



Since these matrices satisfy the two conditions formulated in Theorem the op- 
erator T(A, B) defined on Q is self-adjoint. One can identify the metric graph Q 
with the real line, and under this identification the operator T(A, B) corresponds to 
the operator — -4- sgn(x) 4- defined in L? (R) on its natural domain. 

From the self-adjointness of T(A, B) one deduces that the eigenvalue equation 
can have square integrable solutions only for real spectral parameter. The self- 
adjoint operator T(A, B) can be interpreted as a realization of a system of the 
type 

"-A " 
+A 

which is given with respect to the decomposition £Ul = (£+Ul+) U (£- Ul_). 
The operator T(A, B) generates a unitary group U(t) = e - %tT ( A > B ) which provides 
solutions for the initial value problem 

(i Jj — T(A, B)) u(x, t) = 0, teR, 
u(-,0) = Mo- 
Compared to the time-dependent behaviour of groups generated by (positive) self- 
adjoint Laplacians — A, the group generated by T(A, B) describes a time-reversed 
behaviour on the negative edges and the forward-directed behaviour on the positive 
edges. 

4. Self-adjointness in Krein spaces 

One can consider the same graph described in Example 13.51 but one imposes 
boundary conditions defined by 



.4 



-1 




and B 



These define the non-self-adjoint operator T(A, B) and the self-adjoint operator 
—A(A, B), compare for example [19] . Here —A(A, B) denotes the Laplace operator 
defined in (0J with domain Dom(A(A, B)) = {iJj e V \ Af_ + Bitf = 0}. The 
operator T(A, B) corresponds, again by identifying this graph with the real line, 
to the non-self-adjoint operator — sgn(x)^- with domain H 2 (K) C L 2 (R). where 
H 2 (R) denotes the Sobolev space of order two. A generalization of — sgn(x) from 
intervals to finite metric graphs can be done on the same lines as for — ^ sgn(ir) 4- . 
Define the Krein space 

R = L 2 (Ij), with indefinite inner product [•, — (•, J n ,m') J Hi 

j&SUX 
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where with a slight abuse of notation one sets 

The multiplication with J n)TO is the fundamental symmetry of the Krein space R. 
Since 

J n , m T max = -A and J„, m T min = -A 
holds, where A is the restriction of A to 2?°, one obtains 

Proposition 4.1. The extensions o/T mm in % with J n>m T(A, B) self-adjoint are 
defined exactly by those A, B for which — A(A, B) is self-adjoint in %. 

As already mentioned the self-adjoint realizations —A(A, B) of the maximal 
Laplacian — A are characterized completely, see for example |19| . The eigenvalue 
problem for operators of the type — sgn(x)^- on compact finite metric graphs, 
which have been described here in Proposition 14.11 is considered in [7] . 

Remark 4.2. Analogue to the case treated in Proposition ]^. 1\ generalizations of the 
operator 

-sgn{ X )j-sgn(x)±- 

from intervals to metric graphs can be given. One can search for those extensions 
— A(A, B) of the minimal Laplacian —A with — J n . m A(A,B) self-adjoint in % (or 
equivalently — A(A, B) Krein-space-self-adjoint in Analogue to Proposition ]^. 1\ 
one obtains that the operator — J n ^ n A{A 1 B) is self-adjoint in % if and only if 
T(A,B) is self-adjoint in T-L. Consequently these extensions can be parametrized 
using Theorem \3.1\ 

5. Indefinite form methods 

An approach that has turned out to be very fruitful for sign-definite operators 
in Hilbert spaces is the one using quadratic forms. Representation theorems give a 
one to one correspondence between closed semi-bounded forms and semi-bounded 
self-adjoint operators, see for example (T7l Chapter VI]. The quadratic form defined 
by — -4- sgn(a^) 4± is symmetric, but not semi-bounded. However there are gener- 
alizations of the first representation theorem to indefinite quadratic forms. One 
formulation is given in [12J, for further information on indefinite quadratic forms 
see also the references given therein, in particular the works [10J and [24] . The 
application of results about indefinite quadratic forms to certain indefinite second 
order differential operators on bounded domains is going to be discussed in |16| . 
An analogue construction can be given for certain compact metric graphs. Here the 
result is stated and compared to the previously obtained results, where methods 
from extension theory have been used. 

Assume that (Q, a) is a compact star graph that is one has finitely many internal 
edges X and the initial vertices of all edges are unified in one vertex and all terminal 
vertices are vertices of degree one. This means d-(i p ) = d-(i q ) for i p ,i q € 1 and 
<9+(i p ) d+(i q ), for i p i q . On the compact star graph ((?,§_) one considers the 
spaces H 1 {G) C W and H^\g) C W, where 

H\g) :={i/»eW| Vi(0) = V'i(O), h j € 1} 

and 

Hl(Q) := e H\Q) | = 0, J G l) . 

Define now the gradient operator in % by 

D:Hl{Q)^U, ^V', 
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whose adjoint operator is 

D* : H 1 (G) -> H, ip^-ip'. 

The operator D and the space RanD C H are closed. Since Rani? is closed, it is 
itself a Hilbert space with the Hilbert space structure inherited from H. Hence 



Q: RanD, 



u € Rani?, 
u _L RanD 



is a partial isometry. The adjoint Q* is the embedding from RanD to W. Note that 
(RanD) 1 ^ is the space of constant functions and therefore the operator Q becomes 
more precisely 

Qu = u- — [ u, 

since Q maps onto the orthogonal complement of the constant functions. 

The operators D and D* were used implicitly in [23J to apply form methods to 
Laplace operators on finite metric graphs. The following theorem is a graph version 
of the result that is going to be presented in [16]. 

Theorem 5.1. Let (£/,a) be a compact star graph and 

(a) A e L°°{g;R), that is A, € L°° (Ij ; R) , for each j € I, be such that 

(b) the operator QMaQ* ■ RanD — > RanD is boundedly invertible, where Ma 
is the multiplication operator 

M A : U^U, 4>^A(-)4>. 

Then 

(i) there exists a unique self-adjoint operator C with Dom(C) C Hq(Q) such 
that 

(<P,£il>)H = {</,A(-)il/)n 
holds for all ip 6 Hq(Q) and all ip 6 Dom(C), where 

Dom(c) = e Hl{g) I MaDtP e H\g)}. 

For any ip £ Dom(C) one has Cip — D*MaD%P, the domain Dom(C) is a 
core for the gradient operator D; 

(ii) the operator C is boundedly invertible and its inverse CT X is compact. In 
particular, the spectrum of C is purely discrete. 

Roughly speaking, Theorem 15.11 states that the operator C is a self-adjoint re- 
alization of the formal differential expression — -^A(-)-^ with Dirichlet boundary 
conditions. The proof of Theorem l5 . 1 I relies on the representation theorem for indef- 
inite quadratic forms [T21 Theorem 2.3 and Lemma 2.2]. It is completely analogue 
to the one that is going to be given in [TB] and omitted here. The condition (b) of 
Theorem 15. II can be verified by 

Proposition 5.2. Let (<?,a) be a compact star graph and A(-), ^4(-) _1 G L°°(Q;WL). 
Then QAQ* is boundedly invertible if and only if 

The proof is based on the fact that the space (RanD)- 1 = KerD* is the one 
dimensional space of constant functions in H 1 (Q). It is analogue to the one that 
is going to be presented in [TB] and omitted here. Applying Proposition 15.21 to the 
function 



Jn.m- Q ^ C, Xj I V 



■i. iel- 
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where X = I + U I_ one obtains that QMj n m Q* is boundedly invertible if and only 
if 

a% ~ ai ^ °' 

Then by applying Theorem 15.11 one obtains that the form I„. m given by 

defines uniquely the operator <C n , m = D* J n , m D in H with natural domain 
Dom(£„ :m ) C Hq(Q). Note that the operator £„ jm is a self-adjoint extension of 
T mm and therefore there are operators A and B in /C such that £„. m = T(A, _B). 
More precisely, the operator T(A, B) is defined by Dirichlet boundary conditions 
on the vertices of degree one and at the central vertex v = d—(i p ), i p £ X by the 
local boundary conditions that are given by the matrices 



1 -1 
1 






-1 
1 












Bu 





-1 -1 



where in the last row of B v for each edge in I + stands a +1 and for each edge 
in X_ a — 1 . To paraphrase, these local boundary conditions guarantee that func- 
tions are continuous at the central vertex and that the sum of the outward directed 
derivatives evaluated at the positive incident edges equals the sum of the outward 
directed derivatives evaluated at the negative incident edges. These local bound- 
ary conditions arise naturally from the form approach and therefore they are the 
main example for self-adjoint boundary conditions discussed here. Note that these 
boundary conditions are related to the so-called standard or Kirchhoff boundary 
conditions, compare for example [23J. If A S J and define the standard boundary 
conditions at the central vertex then one has A v = A S J and B v = BffJ n ^ m . 

Remark 5.3. The form approach admits straight forward generalizations to com- 
pact graphs, as long as the operator D: Hq(Q) — > % is boundedly invertible. This 
is equivalent to the invertibility of the Dirichlet Laplacian -A^ = D*D. Hence 
the form approach applies whenever there is at least one vertex of degree one. The 
obstacle that the form approach cannot be used for arbitrary not necessarily compact 
finite metric graphs seems to be due to technical difficulties. 



6. Extension theory background 

One can make one step back and have a look at the problem from the more 
general viewpoint of extension theory. Extension theory deals with self-adjoint 
extensions of closed symmetric operators in Hilbert spaces. Actually it has been the 
starting point of this work to observe that the operator T mm is a closed symmetric 
operator with equal deficiency indices. 

6.1. Classical extension theory. The classical theory developed by von Neu- 
mann describes self-adjoint extensions of closed symmetric operators in terms of 
unitary mappings between the deficiency spaces, see for example Instead of 
the concrete operator J 1111111 consider now two arbitrary closed symmetric operators 
A + and A- in Hilbert spaces H+ and respectively, each with equal deficiency 
indices 



d+(A + ) = d-(A + ) (< ex)), = (< «)), 
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where d±(A±) = dimKer (A* =p i). Then one can study two operators 

T := A + © -A- and A:=A + ®A_. 
Both are closed operators in H = H+ with equal deficiency indices 

d±{A) = d±(A+) + d T (AJ) and d±(A) = d±(A+) + d±(A_). 
In von Neumann's theory the deficiency spaces 

N ± (z) = Kev{A* ± - z), zgC\R 
are important objects. One considers Nt{z) := Ker(T* — z). By construction 



A* + - z 






A*_ + z 



N T (z) = Kcr 

and therefore N T (z) = N + (z) ® N-(—z), whereas for N A (z) = Ker(A* 

N A (z) = Ker 



A* + -z 
A* - z 



holds and therefore N A (z) = N + (z) N-(z). The self-adjoint extensions of T are 
in a one-to-one correspondence to the unitary mappings 

U'\ N T (i) -> N T (-i). 

Analogous to the above the self-adjoint extensions of A are in a one-to-one corre- 
spondence to the unitary mappings 

U: N A (i) -> N A (-i). 

One observes that N-(i) and AT_(— i) are unitarily equivalent, because A_ is a 
closed and symmetric operator with equal deficiency indices. Let J: N-(i) —¥ 
N- (—i) be such a unitary equivalence. The relation between unitary mappings U 
and U' is summarized in the following diagram 



N T {i) = N + (i)(5N^(i) 
u 

N A {i) = N + (-i) © N-(-i) 

and one reads from this that the relation is a one-to-one correspondence, this means 
each unitary map U defines uniquely a unitary map U' and vice versa. This in turn 
gives a one-to-one correspondence between the self-adjoint extensions of A and the 
self-adjoint extensions of T. 



> N+(i)(£N-{-i) =N T {i) 

u' 

> JV+(-i)©iV_(i) = N A (i) 



6.2. Spaces of boundary values. In the context of differential operators the de- 
scription of self-adjoint extensions in terms of boundary values can be more practi- 
cal. Following A. N. Kocubei, see [T51 Theorem 3], for each symmetric operator X, 
defined in a Hilbertspace W, with equal deficiency indices, d+(X) — d—(X) < oo 
there exists a Hilbert space K, of dimension d = d + (X) and linear transformations 
T 1 , T 2 : Dom(A*) -> JC with the properties, 

(1) for any ip, ip E Dom(A*) 

- {v,x*^) n = (rV,r 2 v) K - (rV,rV>/c, 

(2) for any K\,K2 E K there is a ip E Dom(A*) such that T 1 ^ = K\ and 
T 2 ip — K2, and 

(3) ifcpe Dom(A), then T 1 ? = T 2 <p = 0. 
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This can be interpreted as a generalization of Green's formula or simply of integra- 
tion by parts. The triple (/CjT 1 , T 2 ) is called space of boundary values of X. The 
self- adjoint extensions of X are given in terms of unitary mappings U in ]C, see |18l 
Theorem 2 and 4]. More precisely one has that all self-adjoint extensions X of X 
are restrictions of X*, which can be parametrized in terms of U, that is one has 
X = X\j. where Xuip — X*ip and 

DompQy) = {ip G DompT ) | (U - l)rV + i(U + l)T 2 ip = 0}. 

Consider the same situation as above with closed and symmetric operators A± 
each with equal deficiency indices. Let (/C±, T±, be the spaces of boundary 
values of A±. Then the space of boundary values for A is given by (/C,r^,r^) 
with /C = K.+ ffi K- . 



r 1 

1 A 









ri 



and 



r 2 

1 A 









r 2 _ 



which are given with respect to the decomposition K, = K.+ ® K,-. For T the space 
of boundary values is given then by (/C,r^,r|n), again with /C =/C+®/C_, but now 
with 



r 1 - 
i T — 



o 





ri 



and 



r 2 _ 

It - 



r 2 + 
o 



o 

-r 2 _ 



As in the parametrization of von Neumann there is a one-to-one correspondence 
of the self-adjoint extensions of A and T. Given a unitary map U in K, then Ajj is 
the restriction of A* to 

Dom(A c/ ) = {ipe Dom(A*) | (U - + i(U + l)T 2 A tp = 0}. 

With the same U one obtains the extension of T as a restriction of T* with domain 

Dom(T[/) = {ip G Dom(T*) | (U - l)T^tp + i(U + t)T%,ip 
where Dom(T*) = Dom(A*), compare Remark I3.3I 



0}, 



6.3. A radially symmetric example. To come full circle one goes back to the 
starting point - the differential operator given in {1} . Consider the formal differen- 
tial operator r in L 2 (R 2 ) defined by 



-div.A(-) grad it, 



A(x) 



x g 

x G n_, 



where 



0+ = {x G R 2 | ||x|| < Ri} U {x G R 2 | ||a;|| > i? 2 } and 

= {x G R 2 | Ri < \\x\\ < R 2 } 

with R2 > R\ > 0. This ring geometry has been studied in [5]. As A(-) is radi- 
ally symmetric one can transform the operator r to polar coordinates. Using the 
divergence and the gradient operator in polar coordinates one obtains 

Tu(r,8) = ljL r A(r)-^u(r,6) + ^A(r)-^-u(r,6), 
r ar ar r z dt) z 

where 

' -1, re (0,i?i]U[i? 3 ,oo), 
v -1, re (i? 2 ,i? 3 ). 
Using a separation of variables one obtains that r is unitarily equivalent to 



A(r) = 
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where r m are operators in £ 2 ((0, oo), dr) defined by 

, , dr dr , . ., .m 2 -4 _1 , 

T m u m {r) = ~—A(r) — u m (r) + A(r) u m ir). 

or or r 

Denote by ACi oc the set of locally absolutely continuous functions on the interval 
(0,oo). The natural (maximal) domain for r m in £ 2 ((0, 00), dr) is 

Dom(T m ) = {u m G L 2 (0,oo) \ u m ,A{-)u' m € ACi oc anA L 2 ((0, TO ),dr)} 

for m 7^ and for m — one imposes additionally the assumption (or rather 
boundary condition) 

lim[Vrln(r)] _1 uoW = 0. 

9 >Q 

From the general considerations in this section or from the Sturm-Liouville theory, 
see for example [221 122] > one can deduce that the operators 

T m = (r m ,Dom(r m )) 

are self-adjoint. Transforming back one obtains that T = (r, Dom(T)) on the 
natural domain 

Dom(T) = {ue L 2 (R 2 ) \ u G H^R 2 ), -div A(-) gradw G L 2 (n)} 

is self-adjoint. Observe that Dom(T) is a subset of iJ 1 (R 2 ) and elements u G 
Dom(T) satisfy the matching condition for the normal derivatives at the threshold 
of f2+ and fi_ 

lim — u(r + e,8)~ — Um — it(r + e, 6), 
e-s-o+ dr e->o- dr 

where (r,9) G 917+ = <9fi_. 



7. Eigenvalues, resonances and resolvents 

The study of the spectral resolution of the self-adjoint operator T(A, B) is based 
on finding solutions of (t — k 2 )u(-, k) = that satisfy the boundary conditions along 
with certain integrability conditions. Considering each edge separately, one recalls 
that a fundamental system of the equation 

(6) -u"(x,k) = k 2 u(x,k) with k ^ 0, 
u: R — > C is given by e lfea: and e- lfe:E and consequently of 

(7) k) = k 2 u(x, k), for k^O, 

by e _K;r and e K2: . As the two fundamental systems have different integrability 
properties in certain regions of C, one defines the two open quadrants 

Q ={k G C I Re (fc) > and Im (fc) > 0} 

and 

V ={k G C I Re (fc) < and Im (fc) > 0}. 

For fc, k G <2 the functions e lfex and e~ K2: are elements of L 2 ([0, 00)), whereas e~ tkx 
and e KX are not in L 2 ([0,oo)). For fc, k, G the functions e lfea; and e Ka: are square 
integrable. 



SIGN-INDEFINITE OPERATORS ON METRIC GRAPHS 



15 



7.1. Resonances and non zero eigenvalues. A general Ansatz for a solution 
that satisfies simultaneously equation (O on the positive edges and equation on 
the negative edges is the function ip defined by 



(8) 



ip(x, k, in) 



Sj(k)e 

aj(k)e ikx i +I3j(k)e 

a (iK)e~ KX i + f3 3 (in)e KX J 
with k, k £ C \ {0}. With this notation one has 

(r - k 2 )il){x,k,±ik) = and (-A - k 2 )ip(x, k, ±k) 



J G £+, 

j G X+, 
3 G £-, 



0. 



Assuming that £_ 7^ and £+ ^ one has for k £ Q, that ifj(-,k,ik) is square 
integrable. Analogue for k £ "P, tp(-,k,—ik) is square integrable on the external 
edges, as now e fea: S L 2 ([0,oo)). Introduce for brevity the notations 



X+(k) 



'{sj(k)} jeS+ 
{a {k)} jeI+ 



{Sj(iK)} jeS - 

{a 3 (iK)} ]eI _ 



x(k,iK) 



X+(k) 



For the boundary values of the Ansatz function ?/)(•, fc, jk) one obtains the formulae 

/c, in) 

with 



X n ^ m (k,in,&)x(k,ik), if)(-,k,iK)' = Y n , m (k,in,&)x(k,ik) 



X n ,m{k,iK,&) — 



X n (k, a^) 


i^j (fc, a^) 






A m (i/v, a_ 


Y m {in, a_) 



and 



matrices in /C, which are written with respect to the decomposition K, = K,+ 1 
This uses the notation 



X;(fc,b) 



and yi(fc, b) = ifc 







where for Z = n one has to plug in b = a + and for I — m one inserts b = a_ . Actually 
these are the matrices known from the spectral theory of Laplace operators on 
(G+t S+) an( i see for example (2TJ Section 3]. The matrices X n (k,& + ) and 

V n (fc, a,) define operators in K, + and X m (k, a_) and Y m (k, a_) in /C_, respectively. 
The matrix e lfc - denotes the diagonal matrix with entries {e — Sk,ie lkb ' , 

where bi £ b and 5k,i denotes the Kronecker delta. 

Note that there is a function ip(-,k,in) of the form © satisfying the boundary 
conditions Aip(-,k,in) + Bip{-,k,iK) = 0, if and only if there exist coefficients 
x{k,in) £ JC\{0} such that 

(9) Z nim (A,B, k,iK,a)x(k,iK) = 

holds with 

Z n ,m {A, B, fc, ift, a) = AX n m (k, %k 1 a) + BY nm (/c, jk, a). 
Since /C is a finite dimensional space the condition 

det Z ntm (A, B, k, in, a) = 
is equivalent to the condition given in equation (0). 
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Remark 7.1. Let T(A, B) be self-adjoint. Then there are no functions ip(-, k, ik) of 
the form ©, forkeQ with (r - k 2 )tp(-, k,ik) = and ip(-,k,ik) G Dom(T(A, B)). 
Similarly for k G V , there is no such ip(-, k, —ik) with (r — k 2 )ip(-, fc, —ik) = that 
satisfies ip(-,k, —ik) G Dom(T(A, B)). This is due to the fact that tp(-,k,ik) and 
tp(-,k,—ik), respectively would be eigenf unctions to the complex eigenvalue k 2 G 
C\R, which would contradict the self-adjointness ofT(A,B). The question where 
the resolvent R{k) — (T(A,B) — fc 2 ) -1 as function in k admits a meromorphic 
continuation to the real line is related to equation @ . The poles of the meromorphic 
continuation of R(-) are linked to certain singular points of Z n ^ m (A, B, ■, •, a). 

For star graphs, that is for graphs with 1 = 0, the matrix Z ntm (A,B,k,iK,&) 
simplifies to 



Instead of the matrices A and B one can consider the equivalent parametrization 
according to Corollary 13.41 There is an orthogonal projector P and a Hermitian 
matrix L acting in KerP such that with A' = L + P and B' = F J n ,m on e has 
T(A',B') = T(A,B). Denote by lj the eigenvalues of L and by Pj the orthogonal 
projector on the corresponding eigenspace. Then one has for k, k G C \ {0} the 
representation 



Analogue to [T9l Theorem 3.1] one proves now 

Lemma 7.2. Let T(A, B) be self-adjoint. Then zeros of the function 

k h> det Z n>m (A, B, k, ik,a) 

are discrete and there are no zeros in Q. For I. — there are \£\ zeros in Q at the 
most. The analogue statements hold for the zeros of km- det Z„ !m (A, B, in, n,a). 

Proof. The function k i— > det Z n>m (A, B, k, ik, a) is an entire function. This fol- 
lows directly from the definition. Therefore it can be represented as a power 
series. If the zeros had not been discrete, then there would be a convergent se- 
quence {fcjjjgN such that det Z r ^ m (A, B,kj,ikj, a) = 0. This would imply that 
det Z ntm (A, B, k, ik, a) = 0. Therefore it is sufficient to find one point fco G 
C for which det Z n ^ m (A, B, fco, iko, a) ^ 0. If there had been such a fco with 
det Z n ^ m (A, B, fco, ifco, a) = from the open quadrant Q, this would imply the 
integrability of the Ansatz function ip(x, fco, iko). According to Remark 17.11 then 
fco G C \ R would be an eigenvalue of T(A,B), which is a contradiction to the 
self-adjointness of T(A, B). □ 

Putting the pieces together one obtains 

Proposition 7.3. Let T{A, B) be self-adjoint. Then k 2 > is an eigenvalue of 
T(A, B) if and only if there is a coefficient x{k, ik) 7^ with Sj(k) = for j G £+ 
and for the positive square roots fc > of k 2 such that 



The number —k 2 < is a negative eigenvalue of T(A, B) if and only if for the 
positive square root k > of k 2 there is a coefficient x(jk, k) ^ with s 3 (k) = 
for j G £- such that 



Z n ^ m (A, B, k,in) = A + B 



ik 









3 



Z n ^ m (A,B,k,ik,a)x(k,ik) = 0. 



Z n ,m(A, B,in, K,a)x(in, n) — with Sj(re) = 0. 
In particular for 5 = one obtains 
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Corollary 7.4. Let T{A 1 B) be self-adjoint and assume that £ = 0. Then the 
positive square roots k > of positive eigenvalues k 2 are exactly the solutions of 
the secular equation 

det Z n ^ m (A 7 B, k, ik, a) = 

and the positive square roots n > of the absolute values of negative eigenvalues 
— k 2 are exactly the solutions of the secular equation 

det Z niTn (A, B, in, k, a) = 0. 

Remark 7.5. In the presence of external edges only the solutions that are square 
integrable on the external edges are eigenvalues. In general one cannot exclude 
other resonances in the set of singular points of the function k M> Z niTn (A, B, k, ik). 
In all the examples studied the resonances are only eigenvalues, but the author has 
not been able to find a general proof. 

7.2. Eigenvalue zero. The solutions of u"(x) = are the afhne functions and 
only the trivial solution is square integrable on [0,oo). Thus the eigenfunctions of 
T(A, B) to the eigenvalue zero have to be affine on each edge and one gets to the 
Ansatz 



0. 



3 e£+u£- 
■fix,, jex+ux. 



The coefficients appearing on the right hand side are written into the vectors 



and 



(10) 



{a° j } j ex+ 

which are summarized in one vector 



X°_ = 







X° = 



For the boundary values one obtains 
/=X° ro (a) X ° 

with 

^0(a + ) 



4 



and V, ' = r„° m (a)x 



*£(a_) 







which are matrices in JC, that are written with respect to the decomposition /C 
JC+ ® JC-. This uses the notation 



X?(h) = 



and y ; °(b) 



where for I = n one has to plug in b = a + and for I = m one inserts b = a_ . Actually 
these are the matrices known from the spectral theory of Laplace operators on finite 
metric graphs. The matrices X°(a + ) and Y®(& + ) define operators in JC + , and the 
matrices X^(&_) and y^(a_) define operators in /C_. The matrix b denotes the 
diagonal matrix with entries {b}fc,; = 5k,ih. This gives 

Proposition 7.6. Let £ = 0. Then T(A, B) is invertible if and only if 

det(^ iro (a) + J BF„ ro (a))=0. 

The dimension ofKerT(A,B) is equal to the dimension of the space 

Ker(AX® m (&) + BY® m (&)) . In the presence of external edges, zero is an eigenvalue 
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ofT(A, B) if and only if there exists a vector x° with entries of the form (|10p such 
that 

(AXl m (a)+BY r l m (z))x° = 0. 

Remark 7.7. Since the solutions of u"[x) — are affine functions one can use 
the same Ansatz ip° to find solutions of the Laplace equation Au = on the graph. 
Consequently zero is an eigenvalue ofT(A,B) if and only if zero is an eigenvalue of 
the not necessarily self-adjoint Laplace operator — A(A, B). Note that —A(A, B) has 
the same domain as T(A,B), since it is defined by the same boundary conditions. 

7.3. Eigenvalue asymptotics. Assume T(A, B) to be self-adjoint, then its eigen- 
values, according to Lemma 1721 are discrete with only finite multiplicities. Denote 
by < A~i < < . . . the positive eigenvalues, counted with multiplicities and by 
< Af < A^T < . . . the absolute values of the negative eigenvalues, counted with 
multiplicities. Note that in general the eigenvalues of T(A,B) can be embedded 
in the continuous spectrum. One defines the counting function for the positive 
eigenvalues N + (X; T(A, B)) and the counting function for the negative eigenvalues 
N^(X;T(A,B)) by 

N ± (X;T(A,B)) = ]T 1. 

In certain cases a generalized Weyl law holds for the asymptotic behaviour of the 
eigenvalue counting functions. 

Proposition 7.8. 

(1) Assume that (Q, a) is a compact metric graph, that is £ = 0, and let T(A, B) 
be self-adjoint. Then asymptotically the eigenvalue counting functions are 

N±(X;T(A,B)) ~ - ( a M AV2 ' f° rX ^°°- 

(2) Assume that (0,a) is non-compact, but £_ — {£ + = 0) and let T(A,B) 
be self-adjoint. Then the asymptotic formula for iV_(A; T(A, B)) 
{N + (X;T{A,B))) still holds. 

Proof. Suppose that £ = 0. Then it follows from the compact embedding of T> 
H that T{A, B) has compact resolvent and therefore there exists k 2 € R with 
k 2 e p(T(A, B)) n p(T(l,0)). According to the forthcoming Proposition 031 the 
difference 

(T(A, B) - k 2 )- 1 - (T(l, 0) - k 2 )- 1 
is a self-adjoint finite rank operator. To prove the first statement, one applies the 
min — max-principle to the resolvents. As the difference has only finite rank one 
obtains that 

N±(\;T(A,B)) - N±(X;T(1,0)), for A -> oo. 
The operator T(1,0) decouples all edges and defines (positive) Dirichlet Laplacians 
on each positive edge and Dirichlet Laplacians multiplied with minus one on each 
negative edge. The claim follows by summing up the eigenvalue counting functions 
of the Dirichlet Laplacians on the intervals. 

To prove the second statement one takes into account that for a bounded operator 
one can apply the min — max-principle to characterize the discrete spectrum below 
the essential spectrum. For the case of £+ ^ and £_ = the essential spectrum 
of T(A, B) is [0, oo ) and the negative spectrum is discrete. Denote by —/if < the 
largest negative eigenvalue of T(1,0). Then all — k 2 £ (max{- X± , — (X^ },0J are in 
the resolvent set of both T(A, B) and T(l, 0). For such a — k 2 the difference 

(T(A, B) + k 2 )- 1 - (T(l, 0) + k 2 )- 1 
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is a self-adjoint finite rank operator and the negative spectrum of (T(A, B) + fc 2 ) -1 
consists of a sequence of eigenvalues of finite multiplicities that accumulate at zero. 
Therefore one can apply the min — max-principle to characterize the positive eigen- 
values of —(T(A, B) + k 2 )~ 1 as well as the positive eigenvalues of — (T(l, C^ + fc 2 )" 1 . 
As one has again only a perturbation of finite rank the claim follows by the same 
reasoning as above. □ 

For the proof of Proposition 17.81 it is essential to consider only finite rank per- 
turbations. 

Remark 7.9. Let ((?,a) be a compact finite metric graph and assume that T(A,B) 
is self-adjoint and invertible. Then the operator T(A, B) has a compact inverse 
and the eigenvalues of T(A, B)^ 1 can be determined as the successive extrema of 
the Ray leigh- quotient 

(rip, ip) 



ip g Dom(T(A, B)). 



(tV>,tV>) 

This delivers a variational characterizations of the eigenvalues ofT(A,B). 

7.4. Generalized eigenfunctions. An important tool for the study of the abso- 
lutely continuous part of self-adjoint operators T(A, B) and hence for their scat- 
tering theory are generalized eigenfunctions. Assume that £ ^ 0. For I G £+ one 
Ansatz for a generalized eigenfunction is 

'5i 3 e~ lkxi + Slj {k)e lkxi , 
ai 3 {k)e lkx i + f3i 3 (k)e- tkx i 
sij{iK)e^ KX ^ 

^aij(iK)e~ KXj + Pij(in)e KX - 



tpi(x, k, in) 



J e £+, 
j e 1- 



and for I G £. 



(pi(x,k,iK) 



' sij(k) 
a l3 {ky kx i +Pij(k)e 



3 e £+, 

j e i+, 



5ije KX] + sij(in)e KX \ 
y aij(in)e~ KX i + f3ij(in)e K 

with Sij the Kronecker delta and sought after coefficients sij(k),aij(k) and /3ij(k). 
The functions (fi(-, k,in), I G £, are solutions of equation ([6]) on the positive edges 
and of equation (|7J| on the negative edges, but they are not square integrable. For 
notational simplicity these Ansatz functions are written into the \£ \ x \£ UX|-matrix 
<p{x, k, in) with entries 

{<p(x, k, :^ {ipi(x,k,iK)} j , 

that is the (I, j)-th entry of ip(-, k, in) is the restriction of ipi, I G £ to Ij, j £ £ Ul. 
The coefficients for I £ £ + are written into the vectors 



Xl,+ ( k ) = 
and for I G £_ into 

x;,+(fc) = 

Finally one sets 



{<% +s ij (fc)} jef+ 

{ay(fc)} je i + 

{Ai(*)},-er + 

{%•(» }.?'££+ 
{a i O)}j£z + 
{^•(fc)} ie x + . 

Xl{k,in) 



Xl,-i iK ) 



{si J (in)} je£ _ 
{a lj (iK)} jeX _ 



Xl,-( iK ) 



Xl,-( iK ) 



ij(iK)}j 
{otj(iK)} je x_ 



20 



AMRU HUSSEIN 



which is a vector in K, = /C+0/C_. These xi a re summarized into the \£\ x \£ \JI\ 
matrix x(k, in) with entries 

{x(k,iK)} tj ■= {xi(k,iK)} j . 

For the boundary values of the Ansatz functions (pi(-,k,in) defined above one 
obtains the formulae 

ip(-, k, in) = ^L n +m e n,m + X n ^ m (k, ik, a)x(fc, in), 
ip(-,k,iK)' = I n ,m(-k,-iK,)e ntm +Y ntm (k,iK,&)x(k,iK), 



where for brevity 



ikt n 






-nt r 



In 

il„ 



which are related by ifc • 7„, m = I n ,m(k, ik). The notation 







where 



p G {n, m} 



is used, where e„ is a |/Cg , | x |/C|-matrix, e m a |/Cg_ | x |/C|-matrix and consequently 
the block-operator matrix e n<rr , 
conditions 



is a (|/Cg + | + |/Cf_|) x |/C|-matrix. The boundary 
Aipi(-, k, in) + Bipi(-, k, in)' = 0, for I G £ 



are satisfied for all <pi{-, k, in), I G £ if and only if 

{AX n _ m (k, ire, a) + BY n ^ m {k, iK,&))x(k,in) = — (Al„+ m + BI nm {—k, —in))e n ^ m 

holds for the coefficient vector x(fc,i«). If AX n ^ m (k, in, a) + BY nm (k,iK,&) is in- 
vertible one can define the matrix valued transform 

X(k,in) = ~(AX n , m (k,in,&) + BY n , m (k,iK, a)) -1 (Al n+m + BI n ^ n (-k, -in)). 

This is a transform in A and B as well as in k, k. However in the notation the 
dependence on A and B is omitted. When X(k, in) is well-defined one can solve 
the equation for the coefficients x{k, in), which are then uniquely determined by 



(11) 



X(k,in) = X(k,in)e n < 



Note that this way one obtains that x(±/c, ik) is well-defined for all k > 0, except 
for resonances with det Z njm {A, B, ±fc, ik, a) = 0. With the coefficients x(±k,ik), 
where k > such that Z ntTn {A, B, ±fc, ik, a) is regular, the functions 

(fi{-,k,ik) and Lpi{-,—k,ik), I G £+, 

satisfy the boundary conditions At^ + -Be,?;' = 0. By construction they additionally 
solve the equation (r — A: 2 ) ip = 0. Therefore they are denoted as generalized eigen- 
functions of T{A, B) for positive energies. Analogue for I G £- and k > such 
that Z n ^ m (A, B, in, ±re, a) is regular the functions 

ipi(-,in,n) and (pi(-,iK,—K) 

are denoted as generalized eigen] 'unctions of T(A, B) for negative energies. These 
statements are going to be specified, when considering scattering problems related 
to the operator T(A, B) in Section |SJ There the associated wave operators are 
computed in terms of the generalized eigenfunctions of T(A, B). 
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7.5. Coefficient matrix. For a star graph, that is for X — %, the coefficients of 
the generalized eigenfunctions are simplified to become 

X(k,in) = -{A + BI ntm {k,in))~ 1 (A + BI n<m {-k,-in)). 

Given an arbitrary finite metric graph (C?, a) and self-adjoint boundary conditions 
parametrized by matrices A and B one can cut all internal edges in twain and 
substitute each of the pieces by an external edge. Then the graph becomes a union 
of star graphs, on which A and B still define a self-adjoint operator. This motivates 
the definition of the local coefficient matrix of the generalized eigenfunctions £(fc, in) 
for any graph and fixed A,B satisfying the assumptions of Theorem 13. II by 

(12) £(k,m) : = ~( A + BI n>m (k,iK))-\A- BI n>m {k,iK)). 

In short £ is denoted the coefficient matrix. For star graphs £(&, in) — X(h,in) 
holds. The coefficient matrix is an analogue of the local scattering matrix that 
appears in the context of positive Laplacians on metric graphs, see for example 
|19| . and in the case £_ = both objects coincide. With respect to the block 
decomposition induced by £ + Ul+ and £ _ UX_ the coefficient matrix can be written 
as a block-operator matrix 

tfth - \ € ++( k > iK ) £+-(M«) 

where £±±(fc,m) act in IC± and £+±(fc, in) are maps of K,^ to IC±. The poles of 
£(fc,±ifc), counted with multiplicities, are only finitely many, because 
det(A + BI n ^ m (k,±ik)) is a polynomial of degree not greater than n + m and 
therefore it can have n + m zeros at the most. 



Example 7.10. Consider the operator T(A,B) from Examvle \3.5\ Then one ob- 
tains the local coefficient matrix 



£(fc, in) 



1 



fc 2 



k — k -\- 2ink 
2k 2 + link 



2/^ Sz/'l/c 



The local coefficient matrix £ = £(/c, ik) is the k -independent matrix 

_(i + ») -* 

This gives the generalized eigenfunction for positive energy k > 

' e ~ lkx + ie lkx , xe£+, 

(i + i) e - kx , xe£-. 



ipi(x, k, ik) 



Looking at the generalized eigenfunction for positive energy ipi(-, k,ik) one can 

interpret the off-diagonal entry £_| of £ as depth of penetration of the incoming 

wave propagating from the positive edge into the negative part of the graph. The 
diagonal entry C++ describes the scattering behaviour of the incoming waves, in this 

example it is the reflection coefficient. The entries £ |_ and £ have corresponding 

interpretations for negative energies. 

There are some useful symmetries of £(fc, in) summarized in the following 

Lemma 7.11. For complex fc, k ^ and A, B satisfying the assumptions of Theo- 
rem \3.1\ one has 

(1) £(fc,i«;) _1 = £(— fc,— in) and 



(2) £(k,ik) 



1 


0" 




t 0" 





i 




i 



C(-k,iky 



(3) Let be A = P and B = P^Jn^m, where P is an orthogonal projector in K,. 
Then the matrices £(fc,«fc) and £(fc,— ik) are k -independent. 
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Proof, (i) The first statement follows immediately from the definition since 
In,m(-k,-in) = -J n ,mOM«) holds. 

(ii) To prove the second statement, consider instead of A and B the equivalent 
parametrization given in Corollary 13.41 with 



A = L + P and B = P 



1 
-1 



where P is an orthogonal projector and L is a Hermitian operator with L 
P^LP 1 -. This gives 



A - ikP 1 - 
A + ikP 1 - 



1 

-i 

1 ' 

-i 



1 

i 

1 0' 

i 



A + ik 
A-ik 



t 

-i 

t 

-i 



P 2 



ik) = [A + ikP 



For 



and 

<t(k,-ik)* = (^A + ik 
one obtains 

€(k,ik) 



t 
-% 



A - ikP 1 - 



1 

-i 



1 

-i 



P J 



A-ik 



t 

-i 



P A 



1 

i 



<t(k,-ikY 



1 
i 



and therefore using Part (i) of Lemma 17.111 gives the claim, 
(iii) For the proof of the third part of the lemma decompose 



P ± ikP A 



1 

-i 



with respect to the orthogonal spaces RanP and RanP . Denote by 

-l 



P J 



t 

=R 



the inverses of P 



± 



1 

Ti 



P J 



considered as maps in the Hilbert space RanP-K Applying the formula for the 
Schur complement, see for example [33], delivers the fc-independent block-operator 
matrix representation 



€(k,±ik) 



P + ikP A 



1 

=R 



P - ikP A 



1 
Ti 



-IP 





=R 



P J 



-IP 











P J 



P J 



1 

-i 



P 



□ 
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7.6. Resonance equation. The resonance equation © can be rewritten in an 
analogue way to the one known for Laplacians on finite metric graphs, see for 
example [2"T1 Theorem 3.2]. For k, k ^ such that (A + BI n ^ m (k,iK)) is invertible 
the operator Z n ^ m (A, B, k, in, a) can be rewritten as follows 



Z n ,m (A, B , k, in, a) 
--AX n , m (k, in) + BY n ^ m (k, in) 
--(A + BI n<m (k,iK))R+ m {k,iK) + {A 



BI n>m (k, in)) R n m (k, in) 



where 



: (A + BI n . m (k, in)) (1 - €(k, in)T n ^ n (k, in)) . 



\ {^n,m{k, in) + Y n>m [k, in)) . 



become more explicit 



o 



o 

fl+(a_,ire) 



R+(b,k) 



1 
1 







~ik b 



^n,*n(k> S) 



i?„(a + ,fc) 






H-(a_,i«;) 



i? ; -(b,fc) = 




1 

e jfc]i 



and 



T n .m(k,in, a) —R n m (k,in) (i?^ m (/c, iitj) 



T„(a + ,fc) 
T m (a_,m) 



r,(b,fc) = 




e 

e lfe ^ 





ifc b 



Here for / = n one plugs in b = a + and for / = m one inserts b = a_ . Assume that 
k 7^ and that the matrices A + BI n , m (k,ik), A + BI n/m (k, —ik) are invertible. 
Then one obtains the representations 



(13) 



Z nrm (A,B,k,±ik,a) = 

(A + BI n , m (k, ±ik)) (1 - £(fc, ±ik)T n , m (k, ±ik, a)) R+ m (k, ±ik, a) 



7.7. The resolvent. The Green's function for the equation ^— — k 2 ^j u = f, 

f E L 2 (R), is jre ik \ x ~ v '. For computing the Green's function for the problem 
considered here, one has to find appropriate correction terms. These corrections can 

be expressed in terms of the solutions of the homogeneous problem (t-j^z — k 2 ^j u = 
on each edge. The Green's function obtained in this way defines an integral 
operator, which for k 2 E C \ R is the resolvent operator (T(A,B) — k 2 )^ 1 . For 
the computation of the resolvent one follows the guide lines of [2"Tl Section 4] . The 
notion of integral operators is specified in the following definition borrowed from 
there. 

Definition 7.12 ([211 Definition 4.1]). The operator 8. on the Hilbertspace H 
is called integral operator if for all E £UT there are measurable functions 



C with the following properties 



(1) &jji(xj, -)tpj' E L l {Ijr) for almost all Xj E Ij, 
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(2) ip = &p with 

j'esux' 

The (Iu£) x (Iu£) matrix-valued function (x,y) <— > &(x,y) with 

[&( x ,v)]j,j' = &j,j'( x 3,yf) 
is called the integral kernel of the operator A. 

Proposition 7.13. LetT(A,B) be self-adjoint. Then the resolvent 

R(k 2 ) = (T{A,B)-k 2 )-\ for k 2 eC\t, 
is an integral operator with kernel 

r n, m ( x i V, k , + ik ) + rl, m (x, y, k, +ik), k e Q, 

K r n,m( x i y> fc ' ~ ik ) + r l, m ( x > i/> k > - ik )> k^'P, 

where the free Green's function is given by 

o / i • \ \r°(x,y, k) 

n,mv ' — rj^(ar, y, ik) 

where r?(x,y,k), I £ {n, m} are diagonal matrices with entries 

{r?(x,y,k)} p , q =5 p , q e ik ^-^, 



r(x,y,k) 



W n ^ m (k,in), 



W n , m (k,iK) 



W n (k) 
W m (iK) 



where W n {k) = wrl with the identity in K.+ and W m {k) = 4r\ with the identity in 
JC— . The correction term is given by 

r x (x, y, k, in) = -$ n ,m{x, k, iK,)G n>m (k, in, a)$„, m (y, k, iK) T W n , m {k, in), 

where the subscript T denotes the transposed matrix, 

G n ,m,{k, ift, a) = ZK)(-R^ m (/c, &)) Jn,m 

and 

~<5> n (x,k) 
$ m (af,i/c) 



with 



$i(x,k) 



^n,m[Xi k, in) — 

fs t {x,k) 





(p Xl (x,k) tpx t (x,-k)_ 

This notation means that for I = n one plugs in £ n = £ + and T n = 1+ and for I 
•a one inserts £ m = £_ andT rn = Z_. The blocks ipc(x,k), C € {£+,£_,Z+,Z'_} 
are diagonal matrices with entries {(fic(x, k )}i.jec — $l,j{ e 

The proof is postponed to Appendix [21 From the explicit formula in Proposi- 
tion [7TT3] one reads that the resolvent kernel can be defined for all k £ Q as the limit 
of values taken from Q, except for those k £ d Q for which det Z n ^ m (A, B, k, ik, a) = 
0. These exceptional vales are called resonances ofT(A,B) in Q. Analogue for all 
k £ V the resolvent kernel is defined, as the limit of values taken from V , except 
for those k £ dV for which det Z n ^ m (A, B, k, — ik, a) = 0. These are called the 
resonances of T{A, B) in V . This justifies to denote the equations 

det Z njTrl (A, B, k, ±ik, a) = 

the resonance equations for the operator T(A,B). 
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Remark 7.14. From the integrability properties of the integral kernel one reads that 
the resolvent, as a function in k € Q or k € V , admits a meromorphic continuation 
to Q or V. The continuation is possible outside the resonances of T(A,B) in Q 
and outside the resonances ofT(A,B) in V . There the resolvent kernel r(-,-,fc) 
defines an operator 

R £ (k 2 ) : L 2 (G, e ex dx) -> L 2 {G, e~ sx dx), for any e > 0. 

Consequently outside the resonances of T(A, B) the resolvent kernel defines an op- 
erator from L 2 com p(Q , dx) — ¥ L 2 oc (Q,dx), where L 2 com p(Q , dx) denotes the set of 
compactly supported elements of % and L, denotes the locally square integrable 
functions on (0,a). 

Remark 7.15. The values k and in, are square roots of the spectral parameters A = 
k 2 and A = — k 2 , respectively. In the explicit formulae involving k and in it is hidden 
that two different branches of the complex square root are used simultaneously. To 
be specific one uses the branch with Im-yA > and the one with Rc^ > 0. 

8. Scattering 

On the one hand scattering theory gives an interpretation of scattering solutions 
in the context of time-dependent dynamics. On the other hand it gives a description 
of the absolutely continuous part of an operator in terms of perturbation theory. 
As general references for mathematical scattering theory the author recommends 
the books [31| and [2J. 

Here the role of the free operator is played by T(0, 1). This means that in the 
free system all edges are decoupled and on each positive edge one has the positive 
Neumann Laplace operator along with the Neumann Laplace operator multiplied 
by minus one on each negative edge. So, for the self-adjoint operator T(A, B) one 
considers the scattering pair 

(T(A,B),T(0,t)). 

The pre-wave operators W(t) = W{t) {T(A, B), T(0, 1)), where (eE are defined 
by 

W(t) (T(A, B), T(0, 1)) := e ^(^^) e -itT(o,i)_ 
The strong wave operators W± = W± (T(A, B), T(0, 1)) are the strong limits 
W± (T(A, B),T(0, 1)) := a - ( hm^ W{t)P% 0>1) , 

where t \ is the orthogonal projector onto the absolutely continuous part of 
T(0,1). 

Theorem 8.1. Let T(A, B) be self-adjoint. Then the strong wave operators W± ex- 
ist and are complete. The absolutely continuous spectrum o~ ac (T(A, B)) of T(A, B) 
is 

(—00,00), if £ + j^z and £_ 7^ 0, 
[0,oo), if £_ =0 and £+ f 0, 
(-00, 0], if £+ = and £_ f 0, 
0, if £+ = and £_ = 0. 

The multiplicity of {— 00, 0) is \£—\ and the multiplicity o/(0,oo) is \£+\. 

Proof. By Proposition 17.131 the operator 

{T{A,B)-i)- l -{T{Q,\)-i)- 1 

is a finite rank operator and from |31l Theorem 6.5.1] it follows that the strong 
wave operators exist and are complete. □ 



a ac (T(A,B)) 
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8.1. Wave operators. In the stationary scattering theory the wave operators are 
calculated in terms of the resolvents of the perturbed operator T(A, B) and of the 
free operator T(0, 1). The abelian wave operators W| = W£ (T(A, B), T(0, 1)) can 
be computed as 



wis 



lim — ie 



R{\Tie)dE {\)f, 



where dEo(-) is the absolutely continuous part of the spectral measure of the 
free operator and R(X =F ie) = (T(A, B) — (A =p ie)) , see [2, Proposition 13.1.1, 
formula (3)]. Since the strong wave operators exist they agree with the abelian wave 
operators, compare for example |31| . Note that the absolutely continuous part of 
the free operator T(0, 1) is related only to the external edges and can be expressed 
in terms of the cos-transform, see for example |30[ Example 14.8]. The absolutely 
continuous part of the spectral measure of the free operator T(0, 1) is 

dE (X) = ($(dE 0>j (X)). 

For j £ £ + and A > one has 

(dE j(X)f) (x) = cos(VAa; 
and for j £ £_ and A < 

(dEoj(X)f) (x) = cos( v /|A[a; 3 ) / cos( v /|A|y i )/( 2 / i )d % 



^ dX 
] 7rVA Jo 

dX 



■(^%jj)f(yj)dyj 



ttVIA] Jo 

For the cases j £ £+ and A < and j £ £_ and A > one has (dEoj(X)f) (x) = 0. 
Now substitute A > by k 2 — X with k > and A < by — k 2 = X with n > 0. 
Then define on each exterior edge j £ £ the cos-transform by 

2 



and its inverse by 



Set for brevity 



/*(*) = 



fj(k) 



cos{kyj)f{yj)dyj 
cos(kyj)f(k)dk. 



f+(k) = {f j (k)\ , /-(«) = {£(«)} and f(k,n) 



'/+(*)' 
./-(«). ■ 

In the calculation of W± one can interchange the limit e — > 0+ and the in- 
tegration over dEo(X) according to [311 Theorem 4.2.4]. After the substitution 
k 2 = X > one calculates using [311 Definition 2.7.2] 



lim -ieR { k>±ie)dE^)f = lim -ie / r(*,y, >/fc*±fc) 



cos(kx)f + (k) 




ith 



<^ cos(kx)f + (k) > = cos(kxj)fj(k)j 



and denotes the zero on the rest of the components. Analogue after substituting 
—k 2 = X < one computes 



Mm -ieR{-n 2 ± ie)dE {-n 2 )f 



= lim — ie 

e->0+ 



r(a;, y, \/ —k 2 ± ie) 





cos(kx)/_(k) 
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with 

I COs(k2;)/_(k) \ = COs(KXj)fj(K,)j e s_ 

and again denotes the zero on the rest of the components. This calculation 
is carried out in Appendix [B] The result is that the strong wave operators can 
be computed in terms of the generalized eigenfunctions of T(A, B). Note that 
RanP^ Q1) =H £ . 

Proposition 8.2. For f 6 Hs one has 

P OO P oo 

W + f=J2 <Pi(;-k,ik)fi(k)f+Y, V>i{;ik,k)fi(k)f 

,m JO ,,-rn JO 



u leE 

and 



n QQ n QQ 

IGE+ J ° ICE. J ° 

Remark 8.3. From Theorem \8.1\ it follows that the absolutely continuous part of 
T(A,B) is the range ofW± and spectral representations of the absolutely continuous 
part of T(A, B) are given by the transforms 

u± = w±w±. 

8.2. Scattering matrix. As the wave operators exist and since they are complete 
one can define the scattering operator S = S (T(A, B), T(0, 1)) by 

S := W*W_. 

The scattering operator is a unitary operator on the absolutely continuous subspace 
of T(0, 1). From the formulae for the wave operators given in Proposition 18.21 one 
sees that the wave operators decompose into two parts. One part is related to the 
positive absolutely continuous spectrum and the other part is related to the negative 
absolutely continuous spectrum. This observation can be formalized within the 
concept of local wave operators, see for example |31l Chapter 2.2.2]. The local 
wave operators W± (T(A, B), T(0, 1), A) = W± (A) are 

W± (T(A,B),T(0,t),A) := s- lim e «T(A,s) e -itT(o,i) p ac/ A ^ 

t— >±oo 

where Pq c (-) is the spectral projector of the absolutely continuous part of the free 
operator T(0, 1) and A C <x oc (T(0, 1)). For A+ = [0,oo) and A_ = (-oo,0] one 
has, because |A + n A_| =0, that furthermore 

W± = W±(A + ) + W±(A_) . 

The local operators W± (A + ) are 



w+(A+))f= V / <pi(-,-k,ik)Mky- 

IGE+ J ° 

W+(A_))/= ]T / n(;ik,k)fi(k)f 

tr- C JO 



and 



/■OO 

W-(A + )f= V / w(-,k,tk)Mk)f, 

IGE + J ° 

p OO 

W-(A.)f= V / <pi(;ik,-k)fi(k)f 

ir- TP JO 
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For the scattering operator it follows that 

S = W+ (A+)* W- (A+) + W+ (A_)* W- (A_) . 

Let £7° be a given spectral representation for the free operator, then the corre- 
sponding scattering matrix is defined as 

The scattering matrix is unitary. Here one chooses [To to be the diagonal block- 
operator matrix 



U a = 



U+ 

o tv 



where U£f = f + , and U f = /_ 



This is the cos-transform for the absolutely continuous part of the free operator 
T(0, 1), that already appeared in the context of stationary scattering theory. As 
a consequence the scattering matrix is given as a diagonal block operator matrix. 
Using Proposition 18.21 the scattering matrix can be computed in terms of the co- 
efficients of the generalized eigenfunctions. In the following this is carried out in 
detail. 

Write (pi(±k,in,x), I € £. into the \£\ x \£ UX|-matrix &(x, ±k,ik) with entries 
{$(x,±k,ik)} jtl := {ipi{x,±k,in)} j . 
The restrictions of &(x, ±fe, ik) to £ + arc 

{*(*,+M*)W + = Sue-*** +e* kx < X jAk,ik) 

and 

{9(x,-k,ik)} jtlee+ =5i j e^+e- i ^XiA-k,ik). 

Consider the restriction of x(k, ik) to the external edges and denote it by xs,s (k, in), 
which is a \£\ x \£\- matrix with entries 

{x(k, iK)s,s}j,l^S = {x(k, in)} jilee . 

With respect to the division into positive and negative external edges one obtains 
the block structure 



X++(k,iK) x+-(k,in) 
X-+(k,in) x — (M«) 



(14) X s,s(k,in) = 

where Xp,q{k, in) with p, q G {+, — } are the matrices with entries 

{x(.k,iK) Piq } je£p , le s q = {x(k,iK)} je£p l&£q . 

Since the functions e~ tkx and e lkx are on each edge linearly independent there exists 
a \£ + \ x | £ + 1 -matrix C + (k) with entries {C + (k)}ij e £ + = cy(fc) such that 

<f>(-,k,ik)C + (k) = <f>(-,-k,ik) 

holds on the positive external edges. The functions cij{k)tpi(-, k, ik) are generalized 
eigenfunctions as well and therefore their linear combinations 

<Pj(-, —k,ik) — ^""^ c(k)ij(pi(-,k,ik), for x G £ + 

are generalized eigenfunctions, too. Recall that the functions tpi(-, k,ik) and 

—k, ik) are defined uniquely up to a set of measure zero. Hence the relation 

<Pj(;-k,ik)= ^2 c(k)ij(pi(-,k,ik) 
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C{k) 



and hence M c f±{k) = C±(k)f±(k). 



carries over from the external parts to the whole graph (<?,a). Comparing the 
coefficients on the external edges yields 

C+(k)x++(k 7 ik) — Ig, and hence C+(fc)~ 1 = x++{k, ik). 

Analogously one obtains 

(pi(-,ik,k) = ^2 cij(k)(pi(-,ik,-k), 

for an appropriate \£-\ x \£ _|-matrix C-(k) with entries {C_(fc)}j ! ; e g_ = cij(k). 
Considering the restriction of ipi(-,ik,±k) to £ _ yields 

${-,ik,-k)C-(k) = <S>(-,ik,k) 

and hence 

C-(k)x— (ik,k) = l s _ and C_ (fc)" 1 = X — (ik, k). 

This can be used to express W- in terms of W+ as follows, 

(15) W + M c U f = W-f, 

where Mc denoted the multiplication operator with C(-), 

~C+(k) 
C_(fe)_ 

As the absolutely continuous parts of T(l, 0) and T(A, B) are unitarily equivalent 
one can introduce the wave matrices using for example the transform U— from 
Remark 18.31 Here the wave matrices w+ and are defined by 

w+(k)f(k) = U-W+U3f(k) and w.(k)f(k) = U-W-USf(k), 

and they are unitary maps from RanC/o to Ran U- . With the above equation (|15p 
one obtains using UqUoJ = Ufif that 

w* + W-f+ = M c+ f + and w+W- /_ = Af£_ /_ , 

where M£ is the multiplication operator with C±(-). By the definitions of w+ and 
u>_ one has for the scattering matrix 

S(k 2 )}(k) = w* + (k)w-(k)f(k) 

and hence C+(fc) is unitary for almost all k > and as C+(fc) _1 = x++(k,ik) one 
concludes that C+(fc)* = %++(fc,ifc). Analogue one obtains C_(fc)* = X--(ik, k). 
This is summarized in 

Theorem 8.4. Let T(A, B) be self-adjoint. Then the scattering matrix for the pair 
(T(A, B), T(0, 1)) is a \£+\ x -matrix for positive energies A > and a \£— \ x 
|£_ | -matrix for negative energies A < 0, 

S (\) = l X ++^> A>0 ' 

where x++ an d X have been defined in equation (| 14(1 as restrictions of the coef- 
ficient X- 

Recall that the coefficient x has been defined in (fTTj) . Especially one obtains, ap- 
plying Theorem l8.4l to star-graphs, the non-obvious fact that the blocks £ ++ (fc, ik) 
and £ (ik, k) of the coefficient matrix €(k,ik) defined in equation (|12j) . are uni- 
tary for all k > 0, except a finite set. The relevance of Theorem 18.41 arises from the 
fact that it justifies to read the scattering properties of the system directly from 
the coefficients of the generalized eigenfunctions. This is common when considering 
self-adjoint Laplacians on graphs, but it is not self-evident, as the example of a 
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- a 



Figure 1. The graph described in Example 18.5 



Schrodinger operator on the real line with a step potential shows. In this case one 
obtains a relation between the coefficients of the generalized eigenfunctions and the 
scattering matrix as well, but both do not agree in general, compare for example 
the article [13] , which is also helpful for the understanding of the scattering problem 
discussed here. 

Example 8.5. Consider the graph that consists of two external edges £ + = {1,2} 
which are connected by an internal edge I- = {3} of length a = {a}. This means 
that one has the two vertices 9(1) = <9_ (3) and 9(2) = <9+(3). The boundary 
conditions imposed on each vertex are the standard boundary conditions from Ex- 
ample \3.5l They are encoded in the matrices 



A = 



Hence one obtains 



X 2 ,2(k,ik,&) = 



"1 


-1 




















1 








1 





e -ka 









id B = 




1 

oka 





,Y 2 ,2(k,ik,a) 



ik 








-k 


+k 





ke~ ka 


-ke ka 











i k 



and with Zi^A, B, k,ik, a) — AX2,2(k, ik, a) + BY 2t 2{k, ik, a) one obtains the scat- 
tering matrix for k > 



S(k 2 



sn(fc) Sl 2 (fc) 
S 2 l(fc) S 22 (fc) 



itanh(afc) ^^^y 
S^iE^fe) *tanh(afc) 



The absolutely continuous spectrum is [0, oo) with multiplicity two. The negative 
eigenvalues are the zeros of 

det Z2fl{A, B, ik, k,&) — ik 2 cos(afc) 

and there are no embedded eigenvalues. Hence the pure point spectrum 

/((2rn- 1)tt) 2 



a pp {T{A,B)) = {- 



2n 



m G 



consists of infinitely many negative eigenvalues with multiplicity one which accu- 
mulate at — oo. The kernel is zero and there are no further resonances. 



Example 8.6. For the situation considered in Examvle \7.1U\ one reads from The- 
orem that the scattering matrix is 

' +i, A>0, 
-L A<0. 



5(A) 



The absolutely continuous spectrum is the whole real line and there are no eigen- 
values. 
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Example 8.7. Consider a star graph with three edges, E- 
matched together by the boundary conditions 



A 



A 



-1 




and B 



= {1,2} and£_ = {3} 






-1 



A > 0, 
A < 0. 



Note that these local boundary conditions are the ones from Example \3.5\ The 
operator T '(A, B) has absolutely continuous spectrum (— oo,oo), where (0,oo) has 
multiplicity two and (— oo,0) has multiplicity one. Since 

det Z2,i(A, B, k, in) — k + 2ik, 

there are no eigenvalues and no resonances. The local coefficient matrix 
€(k,ik) = £, 

-l/5+(2/5)i 4/5+(2/5)i 2/5-(4/5)i 
4/5 + (2/5)i -l/5+(2/5)i 2/5 - (4/5)i 
4/5+(2/5)i 4/5+(2/5)i -3/5 - (4/5)i 

is k-independent and the scattering matrix is given by the corresponding blocks of 
£ = Xe,£, 

f [—1/5 + (2/5)* 4/5+(2/5)z 
5(A) = \ 4/5 + (2/5)i -l/5+(2/5)i 
1-3/5- (4/5)i, 

Looking at the generalized eigenlunctions and taking into account Theorem 18.41 
one can try an interpretation of the coefficient matrix <£{k,ik). Let Q be a star 
graph and k > 0, then 

• £ ++ [k, ik) is the scattering matrix, that is the entries describe the transmit- 
ted and reflected parts of an incoming wave (scattering matrix for positive 
energies) , 

• €^ (k, ik) can be interpreted as depth of indentation of an incoming wave 

with positive energy k 2 into the negative part of the graph. 

For negative energies ~k 2 < the matrices £__(z/c, —k) and £ ^(ik,—k) admit 

analogue interpretations. 

8.3. Time— dependent Problems. A difficult topic is to consider time-dependent 
problems involving indefinite operators. Thinking of applications to solid-state 
physics, one considers the Schrodinger equation. In the effective mass approxima- 
tion, the effective mass tensor can be negative, too. For any self-adjoint operator 
one can now give the solutions of the initial value problem 

(i§- t -T(A,B))u(x,t)=0, 

u(;0) = U 

in terms of a unitary group 

u = U(t)u , where U(t) = e - itT (A,B) 

is acting in %. The fact that the spectrum is not semi-bounded is no obstacle for 
the construction of this group, similar to the situation, when considering the Dirac 
equation. 

Considering the wave equation or the diffusion equation changes this feature 
completely The lower bound on the spectrum of a semi-bounded operator can 
be interpreted as a measure for the stability of the system. Problems that arise 
when the spectrum is neither bounded from below nor from above can be avoided 
by projecting away the critical parts of the spectrum. Assume that £_ = and 
consider a self-adjoint operator T(A, B) on such a graph. Then T(A, B) has only 
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positive absolutely continuous spectrum. Denote by T ac {A, B) the restriction of 
T(A, B) onto its absolutely continuous subspace. For the wave equation according 
to d Chapter 10.3] with B = T ac (A, B) 1 / 2 the solution of 

d 2 



dt 2 

is given by the group 



+ T ac (A,B) u(x,t)=0, 



u(0) = w , 



d_ 

m 



U = Vq 



'«(*)" 




u'(t)_ 





B -1 shi(tB) 
cos{tB) 





u 







cos(Bt) 
-Bsm(iB) 

acting in an appropriate Hilbert space. All entries are known in terms of the spectral 
theorem using the spectral transform of T ac (A, B), which in turn is explicitly given 
in terms of generalized eigenfunctions, compare Remark 18.31 and Proposition 18.21 
Alternatively one can consider a piecewise defined wave equation 



d 2 


- 


at 7 




a 2 


- -C 


ot 7 


da;- 1 



u(x, t) = 0, 
u(x, t) = 0, 



x e g+, 
x 6 G- 



with appropriate initial conditions. This would correspond to a completely reversed 
dynamics on the negative part compared to the positive one. Such problems - 
although very interesting - are beyond the scope of the approach presented here. 

8.4. Gluing formula for the scattering matrix. Let be given two connected 
metric graphs. Then one can construct a single connected metric graph by gluing 
the two original graphs together along external edges. Each two external edges 
that are glued together in pairs, become an internal edge of a certain length. To be 
more precise let be given two metric graphs {Gi, a^), where Gi = (Vi,X±, £\, d\) and 
(C/2, a 2 ), where Q 2 = (^2,^2) £2,82)- Furthermore let be given two subsets £ 1 C £ \ 
and £2 C £2 of their external edges with \£ x | = \£z\i ancl a bijective identification 
map 

G: £~i^£ 2 . 

Then define a new graph 

G ■= Gi °g G2 

is the graph Q = (V,X,£,d) with V = V x U V 2 and £ = (£i\£i) U ^2^2)- 
The internal edges are X = T\ \JT% UIq, where Xq are the new edges connecting 
the two graphs. One has \Ig\ = \£i\- Any element e g £\ becomes an element of 
i — i e G Tq and one sets d(i e ) — (d(e),d(G(e))). Assigning to each new internal 
length edge the length a,, which is written into a G = {ai}i G i G and keeping the 
edge lengths of X\ and X 2 one obtains a new metric graph ((?,a). 

Let there be furthermore two self-adjoint Laplace operators defined on each of 
the metric graphs (Gj j = 1,2. Then their boundary conditions define naturally 
a self-adjoint operator on the new graph G = Gi°gG2, constructed by gluing 
together the two original ones. One can consider first the scattering matrices of the 
operators defined on each of the components separately. Then the scattering matrix 
of the operator defined on the new graph can be computed in terms of the original 
data of the scattering matrices and the data used in the gluing construction. This 
iterative process is known as the star product. The history of gluing formulas for 
scattering matrices of Laplacians on networks goes back to Redhefier, see [551 US]- 
The star product has been generalized by Kostrykin and Schrader to self-adjoint 
Laplacians on metric graphs, see |19l |2"U] . The considerations presented here are 
based on these last mentioned works. For the history of the star products and 
factorizations of the scattering matrix see also the references quoted therein. Such 
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g 1 



£ l 



G 



I 



-it 



(a) Two disconnetcted star graphs and the identification map G 



Q = Qg 
>-- 

f_ 
a G 



(b) Connected graph after gluing construction with internal edge length ag 

FIGURE 2. Two graphs glued along the negative external edges. 



gluing formulae are only known for one dimensional (singular) spaces. For higher 
dimensions such formulae are not known in general. 

Here a gluing formula is presented for the case of two graphs which have equal 
numbers of negative external edges and all of them are glued together in pairs. 
This allows understanding of the scattering properties of systems, where the neg- 
ative part of the leading coefficient is only compactly supported in terms of the 
generalized star product. 

Consider two finite metric graphs (g 1 ,^, 1 ) and (£7 2 ,a 2 ) with = and 
on each of these graphs the self-adjoint operators T(Ai,Bi) and T(A%, B2), re- 
spectively. In addition let there be a bijective identification G: £_ — > £_. One 
considers now the graph Q — Q\o G Q 2 with Xq C I- and lengths a. L > for 
i € Xq. The boundary conditions {A\,B\) and (^2,-62) act only on V\ or on V2, 
respectively. Therefore one can impose on Q the boundary conditions defined by 



A = 



At 






A 2 



and B = 



Bi 






B, 



Hence T(A, B) is self- 



which are given with respect to the division V = V1UV2. 
adjoint on Q for any choice of a G . 

Let x*(k,ik) be the coefficients of the generalized eigenfunctions of T(Aj, Bj) on 
(Q 3 , a J ), for j = 1,2. With respect to the decomposition into positive and negative 
external edges one writes as in equation (|14l) 

_X-+(k,ik) x--(k,ik) 



X£\ £ i(k,ik) 



and 



t(k,ik) 



x\+(k,ik) 
X-+(k,ik) 



X 2 - 



(fc, ik) 
(fc, ik) 



obtained by gluing together both 
has only positive absolutely con- 



The operator T(A, B) on the metric graph (Q, a) 
along the negative edges, with new lengths a G , 
tinuous spectrum and its scattering matrix is explicitly computable in terms of 
X51 £i{k,ik) and x\^ gs(k,ik), assuming some compatibility properties. To make 
these compatibility assumptions more precise one denotes the critical sets where 
the generalized star product is a priori not defined by 

Si(x\x 2 ,a) :={k>0\ det[l - e - fc * X l-(M*)e~* fi xl-(MA0] = 0} , 
S 2 (x\x 2 ,a) :={k>0 \ det[l - e"*^- (*, ^e-^X-- (k, ik)] = 0} 
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and by 



6(x 1 , x 2 ) = {k > I k singularity of x 1 (fc, ik) or x 2 (k, ik)} 



The generalized star product is denned in the article and is denoted by * p . The 
gluing formula is described in the following 

Proposition 8.8. Let T(A\,B\) and T(^4 2 ,i? 2 ) be self-adjoint operators on the 
graph {Q ,a 1 ) and on (5 2 ,a 2 ), respectively, where p = \£_\ = \£-\- Let fur- 
thermore be T(A,B) and (G,&) as described above. Then for all k > with 
k iX >§) U "2(x >X >&) ^ Q(x 1 7X 2 )} ^ e scattering matrix of the pair 

(T(A, B), T(0, 1)) is given in terms of the generalized star product as 

r„-fci 

^(fc 2 ) = Xe,e(* 5 **) *p ^fe G > k)xl i£ (k, ik), V(& G , k) 



e -ac 




Carrying this out gives 



S(k 2 ) 



S2i(fc) s 22 (fc) 



where the blocks Sy(fc), i,j € {1,2} o/ £/ie scattering matrix are matrices with 
entries 



{ s ij(k)} nm — {(sn,m(fc))}„ e£I+ , mE ^ 



One /ias 



«U =X+ 

.2 



4- 



S21 
S22 



[1- 



,-fci 



[1- 

2 „-fc 



x l_ e -fca Gx 2__ 



]-V_+, 

]-l e -fea Gx 2 



=X+++X+-e fc ^ x i_ e -fca G 



[1- 



,-fc 



a. 



«X--e rcaG X--] X e fc - G X-+, 
is s/iori for x± ±{k, ik), j = 1,2. 



a. 



where the k-dependence is omitted and x± ± 

The proof is completely analogue to [50] and omitted here. It is based on the 
fact that the generalized eigenfunctions of the new problem can be obtained as 
linear combinations of the generalized eigenfunctions of the two original scattering 
problems. The restriction imposed on the energies for which the gluing formula is 
valid is due to the fact that, in contrast to the case of self-adjoint Laplacians on 
graphs, the relation between eigenvalues and resonances has not been clarified yet. 

The formula in Proposition l8.8l shows that for positive energies the matrix V(a,, k) 
damps the scattering waves on Xq down. Let be k n > with 

k n i {S 1 (x 1 ,x 2 ,a)US 2 (x 1 ,x 2 ,a)U6(x 1 ,X 2 )}, neJS, 
and k n — > oo for n — > oo and assume that the limits 
(16) lim x++(k n ) and lim x++(M 

n— too n— >oo 

exist. As a direct consequence of Proposition 18.81 one obtains that then 

lim s 2 \{k n ) = 



lim si 2 (fc„) = 0, 

n— >oo 



and 



lim sn(k n ) = lim x+ + (k n ), lim s 22 (fc„) = Hm x++(M 



which shows that the scattering matrix resembles for large energies the direct sum 
of the scattering matrices of the two building stones. However it is not clear that 
the limits in (JTBJ do always exist. 
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Remark 8.9. Gluing together positive as well as negative edges can be done in two 
steps. First one glues the negative edges together and then one considers this graph 
and an auxiliary graph consisting of two edges with continuity boundary conditions. 
This is used to connect two positive edges. More precisely after having constructed 

S(k 2 ) = XsAK ik) * p V(& G , k)xl £ (k, ik) 



according to Provosition \8.8\ one considers the auxiliary graph Q aux — (V,£,d), 
\V\ = 1, \£ | = 2 and the Laplacian defined on it by the boundary conditions 



-1 




and B 



This operator — A a 



-A(A,B) has the scattering matrix S aux {k 2 





- e ika o" 




"0 


1' 


*2 


1 




1 






Now one glues together two external edges of Q and replaces them by an internal 
edge of length a. The scattering matrix of this new problem is 



S(k 2 



Example 8.10. Consider again Examvle \8.5l The graph described there can be 
obtained also by gluing two star graphs, where each is a copy of the graph given 
in Example \3.5i These are glued together along the negative edges, which gives a 
new negative internal edge of a certain length a = {a}. The coefficient matrices 



X++(k,ik) x+-(k,ik) 
X~L + (fc, ik) x 1 - {k,ik) 



(1 + *) 



(!-<)' 



see Example \ 7.1(A Using the gluing formula from Proposition 
length a one obtains 



for J = 1,2, 

with new 



sn(k) =x 1 ++ +X 1 +-e- ka x 2 - 



-ka 



and 



=i tanh(fca) 



S2i(k) =xl-e- ka 
1 



[i-x 1 



~ ka x 2 - 



-ha~\ 



-ka] 



cosh(ka) 

This is well defined for any k > 0. For symmetry reasons one has S22{k) — Sn(fc) 
and S2i(k) — S\2{k). Applying the gluing formula yields the same result as the 
direct calculation in Examvle \8.5[ Note that the limits 

0" 



lim S(k 2 ) 



and lim S(k 



exist. For small energies there is full transmission, whereas for large energies there 
is full reflection along with a phase shift. In addition, the scattering matrix resem- 
bles for large energies the direct sum of the scattering matrices of the two building 
stones, compare Example I #.61 This is a feature not known for scattering systems 
involving only semi-bounded operators. 

8.5. Scattering and cloaking phenomena? In [5j the light propagation through 
metamaterials is described using regularizations of an indefinite operator 



= — div A(-) grad-i 



A(x) = 



x € f2 + , 
x € 0_, 
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where 

n+ = {x G R 2 | ||se|| < i?i} U {x G R 2 | ||x|| > i? 2 } and 
= {x G R 2 | i?i < < R 2 } 

with i?2 > Ri > 0. The quantity computed there for the measurement of the 
light propagation is the Dirichlet-to-Neumann-map for certain exterior domains. 
In the setting of the Sturm-Liouville theory the Dirichlet-to-Neumann-map is 
given by the Weyl-Titchmarsh m-function, which is also closely related to the 
scattering matrix. A generalization of the Weyl-Titchmarsh m-function to partial 
differential operators on exterior domains has been given recently in As in the 
one dimensional case the generalization of the m-function is again the Dirichlet- 
to-Neumann-map. One can suspect that the generalization of the m-function is 
also closely related to the scattering matrix and to inverse scattering problems. 

The analysis of the Dirichlet-to-Neumann-map in [5] leads G. Bouchitte and 
B. Schweizer to the interpretation that in the considered situation a cloaking phe- 
nomenon is taking place. Cloaking simplified means that there is an object in the 
system that can be made invisible for an observer looking at the system from out- 
side. Of course, one can make any object invisible by putting it into a kind of box 
and decoupling the observer from the object of interest inside. Cloaking means 
rather that the measurement of the observer is the same in both cases, when the 
object is present and when it is not. 

Since at least in dimension d = 1 the Dirichlet-to-Neumann-map which is the 
measured quantity in [5 J is related to scattering properties of the system it is reason- 
able to consider also scattering problems involving indefinite operators. The gluing 
formula given in Proposition 18.81 allows the qualitative analysis of the scattering 
properties of a quasi-one-dimensional system involving compact components, on 
which the symbol of the operator is not sign-definite. In particular one can observe 
that even very small negative components have strong influence on the scattering 
matrix and yield to new features as the one discussed in Example 18. f 01 However 
the author has not been not able to identify cloaking phenomena in the system 
considered. For a further investigation one can think of studying the asymptotic 
behaviour of the scattering matrix for small and for large energies. It is not self- 
evident that these limits exist in all cases. 

9. Summary and outline 

An explicitly solvable model problem for second order differential operators with 
(discontinuous) sign-changing but invertible symbol has been introduced. All self- 
adjoint extensions have been parametrized using methods from extension theory. 
The crucial observation of this note is that there is a one-to-one correspondence 
between the parametrization of the self-adjoint realizations of the model operator 
t which acts as "plus Laplace" or "minus Laplace" and the parametrization of the 
self-adjoint realizations of the positive Laplacian. This can be used to distinguish 
one realization among the set of realizations. This is done in two steps. First 
one considers the natural self-adjoint realization of the Laplace operator. This is 
defined by the so called standard or Kirchhoff boundary conditions and it can be 
parametrized using a unitary map in the deficiency spaces of the minimal Laplacian. 
Taking now the same unitary map to parametrize the extensions of the model 
operator r gives a distinguished self-adjoint extensions of T mm . It turns out that 
this realization is exactly the natural realization of r, which is also obtained using 
indefinite quadratic forms. Natural means that the domain of the operator is the 
natural domain for the composition —-j-A{-)-^-. 
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The conjecture is that this holds in more general situations and that this allows 
to show the self-adjointness of sign-indefinite operators of the type —divA(-) grad 
in their natural domain in L 2 (fl) with A(x) elliptic for x 6 il + and —A(x) elliptic 
for x 6 f2_, where = Q + U f2_. In a very particular radially symmetric situation 
this has been verified. The scheme presented here at least gives a distinguished 
self-adjoint realization of such operators and for the model problem with finite 
deficiency indices the conjecture is indeed true. 

The spectrum of the self-adjoint realizations of the model operator resembles the 
spectrum of —A and +A, where —A denotes the Laplacian on certain subgraphs. 
Using explicit formulae for the resolvents, the wave operators and the scattering 
matrix have been calculated. Even if there are only very small components where 
the symbol of the operator is not sign-definite this has strong influence on the 
qualitative behaviour of the scattering matrix. Particularly the scattering matrix 
exhibits new features and does not just resemble the scattering matrix of a direct 
sum of sign-definite operators. 

It constitutes an open problem whether the sign-indefinite operators that are 
considered in this paper arise naturally as limits of certain operators on domains or 
manifolds of dimension greater than one when the "thickness" shrinks to zero. For 
Laplacians on graphs there are such results, see [51lll| and the references therein. In 
the mentioned articles it has been shown that some self-adjoint Laplacians on finite 
metric graph are obtained as the limit of a Laplacian on a "thick graph" which is a 
thin neighbourhood in R n of the original graph. For Laplace operators on twisted 
tubes similar results have been obtained recently, see [22] . In particular it has been 
shown that the limit of such operators yields certain Schrodinger operators on the 
real line. 



Appendix A. 

Proof of Provosition \7. I'M In order to prove that the kernel r(-,-,k) defines the 
resolvent operator (T(A, B) — fc 2 ) -1 one has to check: 

(i) (T max - fc 2 ) J g r(x, y, k)ip(y) = ip{x) for all ip € H, 

(ii) R(k 2 )tp = J g r(x, y, k)ip(y) € Dom(T(A, B)) for <p e U and 

(iii) r(x, y, —k) — r(y, x, fc)* holds for Im fc > 0. 

The statements (i) and (ii) prove that R(k 2 ) is the left inverse of T(A, B) — fc 2 . To 
prove that it is also the right inverse it is sufficient to verify (iii) . The proofs of (i) 
and (ii) are given for fc 6 Q with (iii) the claim carries over to fc E V. 

Proof of (i) : For ip £ % with ipj 6 (Ij , C) for every j 6 £ UX, one uses that 
the Green's function to the problem is known and one has 

<Pj(yj)dyj =<fj(xj), for fc e Quv, 
(fij(yj)dyj =tpj(xj), for fc e Q and 

Vj{Vi) d yj =l Pj{ x j), for fc e . 

Since the kernel r„ m (x, y, k, in) defines a bounded operator in H, by continuous 
continuation from 

v = cre.c) 

jesux 
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to Ti it follows that for k 6 Q and for all cp € % 



\T r ' 



holds. Observe furthermore that 
(T max - k 2 ) 

which completes the proof. 
Proof of (ii): Observe that 



r°(x, y, k, ik)<p(y)dy = tp(x) 



r 1 (x,y,k,ik)ip(y)dy = 



and hence 



r°(x,y, k,iK)ip(y) € V and / ^fay, k,in)ip{y) € T>, 

Jg 

R(k 2 )p = I r(-,y, k)<p(y) eV = Dom(T max ). 



Let be k £ Q and consider the traces for the free Green's function r°(-, •, k, in). 
One has 

r (x, y, k, in)f{y) =R+ (k, in, a) -1 J n , m f{k, in), 



r (x,y,k 7 iK)f(y)) =I n ^ m {k,in)R+ m (k,iK,&) 1 J n , m f( k ^ K ) 



with 



f(k,in):= / $ n , m (y,k,in)W n , m (k,in)f(y)dy. 



Since the free Green's function decouples the positive from the negative edges the 
above statement follows already from the corresponding calculation for self-adjoint 
Lapalce operators, compare [2"Tl Proof of Lemma 4.2] 

For the correction term r (■, •, fc, in) observe that for an appropriate vector / one 
has the traces 







e lkx e 



' fs~ 




"10 




' fs~ 






1 1 




h- 













e 











ik 







' fs~ 


h- 







ik —ik 













—ike lk - ike~ lk - 







and therefore 

®n, m (x, k, in)f = X n>m {k, m)f, ($ n ,m(x, k, in)f^J =Y n<m (k,iK)f, 

where / is here short for f{k,in) 1 compare also [2"Tl Proof of Lemma 4.2]. This 
gives the traces of the correction term 

/ r 1 {x,y 1 k,in)f{y) = X, hm (k, iK)G„ im (fc, in, a)/ 
±G 

and 

r 1 (x,y 7 k,iK)f(y) I = Y n ^ n (k 7 iK)G n , m (k,in,a,)f. 
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Together one has 



r°(x, y, k, »*)/(») + B {^f r°(x, y, k, 
+aJ r 1 (x 7 y,k 7 iK)f(y) + b(^J r : (x, y, k, ilt)f(y)j 



—AR^ m (k 7 iK,&) 1 J n ,m/ + BI nim (k,itt)Rn^ m (k,iK, a) 1 J n ,mf 
~\~AX n ^ m (k, iK)G n ,m{k, in, a)/ + BY n ^ m (k, in)G n , m (k, in,&)f 

where one has used the definition of G n ^ m (k,iK, a) in Proposition 17. 1^1 It follows 
that 

i?(fc 2 )/ - / r(;y,k)<p(y) e Dom(T(A, B)). 

Proof of (iii): As -^Q lk \ x ~y\ is the Green's function for the problem on i 2 (R), 
the symmetry holds for r°(x, y, k, ik) with k € Q. It remains to prove the symmetry 
of the correction term r (x, y, fc, ifc). Substituting formula (|13p into the formula for 
the correction term for appropriate k, k gives 

r 1 ^, y, fc,i«) =$ n , m (x, fc, iK)(i?+ m ) _1 (fc, i«,a)[l - £(fc, in)T n , m (k, in)]- 1 o 

oX(k, iK)(R+ m y 1 (k, in) J„, m <Iy m (y, fc, iK) T W„ !m (fc, ik). 

Notice that W n , m commutes with Rf L m as well as with $„ im and J„. m . Consider 
now 

r(x, y, k, ik) for k e Q. Using for the diagonal matrices l n ^ m and J n ,m the equalities 

^n.m^m — Ri,m and In,mJn.m — ^nm 

one gets 

{ [1 - £(fc, ik)T nim (k, ifc)]- 1 £(fc, ifc)W„, m (fc, ifc) J„, m }* 
= {[1 - C(A, ik)T n<m (k, ik)]- 1 <L(k, ik)W n , m (k, ik)J n , m }* 
=J n , m W n , m (k, ik)* €(k, ik)*[t - T n>m {k, ik)* €(k, ik)*}- 1 

= —=.I nm J nm ^L{k^ik) [1 T nm (k,ik) ] 
—Ik 

= -|~[l - U:, m £(*, ik)*)I n , m {T n>m {k, iky)]-\ll m £(*, ik)*) 
—Ik 

Using the formulas given in Lemma 17.111 one continues as follows 

-^=[1 - (i:, m £(fc, ik)*)I n , m {T^ m {k, ik)*)]-\ll m €(k, ik)*) 

= -5=[l - £(-*, ik)T n , m {k, ik)*]- 1 £{-k, ik)I* m 
—2k 

= [1 - €(-k, ik)T n ^ n (-k, ik)]- 1 <t(-k, ik)W n , m (-k, ik) J„ >m 
=[1 - £(-*, ik)T, hm (-k, ik)}- 1 £(-k, ik)W n , m {-k, ik)J n , m . 
Furthermore one has 

R^ m (k,ik,&)* — R^ m (—k,ik,a.) and §(x, k,ik)* = $(x, —k,ik), 
compare also [21, Proof of Lemma 4.2]. Putting the pieces together one obtains 
r 1 (fc, ifc, y, x)* = fe, ife, a;, y) for fc € Q 
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which proves the claim. 



□ 



Appendix B. 



Continuation of the proof of Proposition \8.Sl The computation of the wave opera- 
tors is supplemented. One computes for k 2 = A > 



lim ± ie) d E (k) f = hm -ie [ r( X , y, V^TTe) 



cos(fccc)/ + (fc) 




dy 



with {cos(fca;)/_|_(fc)}j g f , = cos(kxj)fj(k)j^s + and with is meant the zero on the 
rest of the components. For — k 2 = A < one computes 



lim -ieR(-K> ± ie) d E ( K )f = lim -ie [ r( X ,y, ^^e) 





cos(nx)f- (k) 



dy 



with {cos(Ka;)/_ = cos(ra j)fj{^)j&8- an d again denotes the zero on the 
rest of the components. 

Fix first the branch of the complex square root 

v 77 : C\ [0,oo) -> C+ withlm^X), 

where C + = {z £ C | Imz > 0}. Consider the limit values 



= \j — k 2 + ie and n t = \/—k 2 — ie. 
Taking the limit e —¥ 0+ for k > and k > 0, respectively gives 



lim fcj~ = k, 



lim fc, = — fc, 



lim kJ — in and lim ft, = in. 

e->0+ c^0+ 

Observe that for e — > 0+ with fixed k > and k > one has asymptotically 



lim (kJ — in) ~ — , 



lim — in) ~ — , 



lim ffcj" — fc) ~ — and lim (k~ + k) 

e->0+ V y 2fc e-s-0+ V ' 



2k 

The following auxiliary calculations are needed, 
* k t\ x -y\ C os{ky)dy = 

/oo 



2fc 



ze / e 
'o 





~ t£ / „ikt{x-y)-iky , 

2 7o 2 

_ 16 f i(k-k+)x 



2i(fc- fc+) 



~ le -ikfx ( _ e i(k++k)x 

2i{kt + fc) V 



i Ze c ikfx ( e -i(k+ + k)x _ A , ie c -ife+x f_ e »(fc+-fc)aA 

-2i(fc+ + fc) V / 2i(fc+ - fc) V / ■ 

In the limit e — > 0+ all summands of the right hand side except the last one vanish. 
The last term becomes —ike~ lkx . Consequently 



lim — i 

e->0+ 



ik+l 



cos(ky)dy = —ike 



— ikx 
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Analogously one obtains 



r-OO 

ie I e lk ^ y cos(ky)dy = 
Jo 

r-OO - r-OO 

_ / e ikfy+iky _^ ~ le I & ik+ y - ik V = 

' Jo 2 J 



" (-D + ^TTl-D. 



2i(kt + k) x ' 2i(k+ - k) 

In the limit e — > 0+ only the second addend of the right hand side remains and 
takes the value — ik. Consequently 

r-OO 

lim — ie / e lk ^ v cos(ky)dy = —ik. 
Jo 

Now the same calculation with k~ yields 

r-OO 

-ie e tk < |x_!/| cos(ky)dy = 
Jo 

■ r-X ■ r-OO 

J e ik~ (x-y)+iky _|_ J & ik~ (y-x)+iky ^ 

- r-X ■ r-OO 

+ J e ik~(x-y)-iky _|_ y gifc^ (y-x)-iky^y _ 

~ l€ c ik7x ( e i(k-k~)x _ A , ~ le c -ik7x (_ e i(k-+k)x\ 

2i{k-K) v j 2i{K + k) v y 



-2i(k7 + k) v y 2i(kr-k) v y 



In the limit all terms on the right hand side except the second term 



-e 



2i(fc e " + fc) 

vanish, which becomes —ike lkx . Consequently 



Analogously one obtains 



r- OO 

lim -ie / e ik < |x_!/| cos(ky)dy = -ike ikx . 
Jo 



/■OO _ _• f- OO _ ■ r-OO 

-ie e lk - y cos(ky)dy = — / e ik * y+iky + — / e ik ^ y ~ iky dy 
Jo 2 J 2 J 



(-1)+ * ,, (-!)■ 



2i(fc e " + fc) 2i(fc e " - fc) 

In the limit e — > 0+ only the second term remains and takes the value —ik. Con- 
sequently 

r-OO 



poo 

lim —ie / e lk ' y cos(ky) = —ik. 
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Taking into account the above auxiliary calculations yields 

dy 



lim — ie / r(x, y, \Jk 2 + ie) 
e— >-0+ 



cos(kx) f (k) 




: lim — ie / r°(x,y, \/ k 2 + ie,i\/k 2 + ie) 

e->0+ 



+ lim —ie / r x (x, y, \/ k 2 + ie, iy k 2 + ie) 

e— ¥0+ 



cos(kx)f(k) 


cos(fcx)/(fc) 




dy 
dy 



-ikx . 





+ 2^n,m{x , k, ik)G n ^ m (k, ik) 



f(k) 




- ipi{x,k,ik)fi{k). 



Analogue 



lim 



e->0+ 

+ lim 

= 1 \e ikx j{k)~\ 1 

2 [ 

=^ X! fi(x,-k,ik)fi(k). 
1&S+ 

Using similar calculation gives 



r(x,y, \/k 2 - 


-ie) 


cos(fcr)/(fc) 



r°(x,y, \/k 2 


- ie, 


— iy fc 2 — ie) 


-r x {x,y, \Jk 2 


- ie, 


— iy k 2 — ie) 



dy 

cos(kx)f(k) 


cos(fcx)/(fc) 




dy 



+ 2 ^n,m( a '; — ik)G n _ m (—k, ik) 



f(k) 




lim — ie / r(x, u, v— k 2 + ie) 

e-S>0+ ' 







and 



lim -ie [ r(x,y,y^Te ) 



cos(kx)/(/«) 





cos(kx)/(k) 



dy = ^ X! fi(x,ik,~k)fi(k) 



Using 



/oo 
Hm + -ie.R(A T ie)dE (\)f, 



one obtains that the kernel of the wave operators is given in terms of the generalized 
eigenfunctions. □ 
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